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Abstract 



Recently, the one dimensional model of N spins with S = \ on a circle, interacting 

with an exchange that falls off with the inverse square of the separation: 

-f^iSE = J2i^j pv ■ h-j Tp(Sj • Sj — \), or ISE- model, has received ample attention. 

L n Sin^ 7TJJ 

This model was introduced by Haldane and Shastry in 1988. Its special features 
include: relatively simple eigenfunctions, non-interacting elementary excitations that 
obey semionic statistics (spinons), and a large "quantum group" symmetry algebra 
called the Yangian. This model is fully integrable, albeit in a slightly different sense 
than the more traditional nearest neighbor exchange (NNE) Heisenberg chain. The 
simplicity of the ISE-model warrants its paradigmatic role as the semionic analog of 
the free boson or fermion gas. 

This thesis comes in two parts. Part I (chapters 1 and 2) deals with the integra- 
bility of the model, and Part II (chapters 3 and 4) discusses some of its applications. 
Chapter 1 introduces the model and presents the construction of a subset of the 
eigenfunctions. The other eigenfunctions are shown to be generated by the action of 
the Yangian symmetry algebra of -f^isE- The Yangian is derived from a transfer ma- 
trix, thus establishing a relation to more common exactly solvable models. Chapter 2 
presents a method to construct the set of constants of the motion of the ISE-model. 
Normally they derive from the trace of transfer matrix, but here they are obtained 
by considering the deformation of the spin model into a dynamical version. 

The ISE-model is tractable enough to obtain its zero-magnetic-field dynamical 
structure factors. Chapter 3 attempts to extend this to a non-zero magnetic field, 
where, due to the presence of spinons in the groundstate, more complicated excitations 
contribute: in low fields just up to two extra spinons, for high fields two oppositely 
moving magnons. We discuss the relation to the more complicated NNE-model. 
Finally chapter 4 illustrates how a recently conjectured new form of Off Diagonal 
Long Range Order in antiferromagnetic spin chains can be reinterpreted as a spinon 
propagator in the ISE-model, and verified numerically. We briefly comment on its 
relevance to stabilizing superconductivity in the layered cuprates. 
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Introduction 



In this thesis, we will confine ourselves almost exclusively to one space dimension, 
and mostly one specific model that hasn't been realized experimentally yet. The ex- 
perimental side of other models that we will discuss is going to receive little attention 
either. Given these serious restrictions, one may wonder to what extent this work 
should belong to the field of Condensed Matter Physics. We will try to address this 
concern in three steps. 

Why one dimensional physics? The fabulous thing about one dimensional Quan- 
tum Mechanics (or equivalently two dimensional statistical mechanics) is that there 
exist interacting many-particle models that are exactly solvable. Unlike in higher di- 
mensions, particles cannot evade each other and have to scatter. For special choices of 
the interaction potential, particles can be made to propagate in plane waves between 
scattering events. During a collision they exchange their plane wave momenta. The 
coefficients of both sets of plane waves (before and after) are related to each other by 
the phase shift, picked up by the two particles that scattered. This heuristic attempt 
at formulating a wavefunction is called the Bethe Ansatz [0, 0]. The Bose gas with 
delta-function interactions is an example where the Ansatz works. 

We use the word solvability in the sense that the corresponding Schrodinger dif- 
ferential/difference equation can be reduced to a system of algebraic equations that 
can be solved, at least numerically 1 . For most models the latter is indeed all one 
can do analytically. This brings us to the second point: why we discuss one specific 
model, viz. the model with N spin-| particles distributed equidistantly on a circle, 

1 E.g. for the Bethe Ansatz, the plane wave momenta are related to each other and to the phase 
shift via a set of nonlinear equations that can be solved numerically for finite N and analytically in 
the thermodynamic limit. 
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coupled with an exchange interaction that decays proportional to the inverse squared 
distance between the spins: 



N 



HlSE = ^2 



^[(f) 2 sin(^ 



S,; • S,-. 



2 kjj yjj 



It is also called the Inverse Square Exchange or ISE-model. In this model solvability 
means a great deal more: not only are a lot of the physically relevant eigenstates 
known analytically (and the others can be generated via the action of ladder opera- 
tors), they are also of a relatively simple form, so one can actually calculate things 
like structure factors in closed form. These are averages of local operators in a certain 
state. To date this is still impossible in the Bethe Ansatz style models. 

The benefits don't stop here; the elementary excitations of this model, so-called 
spinous, are spin-| particles that obey semion statistics, that is to say, for every two 
spinons we create, there is one less orbital available for the next spinon. For fermions 
this would be 2 and for bosons this would be 0. This is an example of the exclusion 
statistics interpretation of fractional statistics [0], as opposed to the usual exchange 
statistics of the two dimensional fractional quantum hall effect (in which case both 
formulations can be shown to be consistent). The long range spin model is not the 
only place where one finds fractional excitations, but contrary to spin models such as 
the traditional Nearest Neighbor Exchange (NNE) Heisenberg model, the ISE-model 
describes a free gas of these spinons. 

In the continuum limit there exists another solvable interacting model, the Suther- 
land model [0], of M bosonic particles on a circle of circumference N, repelling each 
other via a ^ interaction, which has fractional excitations. 

M M i 

#suth = -EC + A (A-i)£' 



fe[(fsin(^r) 



2 



As a matter of fact, one can select their statistics, by varying the coupling constant in 
the potential energy. For a particular value of that coupling A = 2, (corresponding to 
exactly semion statistics) we can set up a mapping to the ISE spin model. The map- 
ping identifies a number of low lying states in the continuum with a number of max- 
imal S z states in the spin model. The correspondence is not trivial in that all other 
(missing) spin states are obtained by the action of some simple symmetry like SU (2) 
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or SU(2) x SU{2) (as in the NNE and ID Hubbard models respectively). Recently 
it was found that the "missing" states are generated by acting on the "continuum" 
wavefunctions with the much larger symmetry algebra called the Yangian, Y^s^), a 
quantum deformation of the SU (2) Lie group. Unlike the case of the Quantum Affine 
Algebra extension of SU(2): U q (sl2), the Yangian is not obtained by deforming the 
Lie group commutation relations, but rather by adding new generators. The Yangian 
algebra is well studied in the mathematics literature and was first introduced into 
physics in quantum field theory [0]. The spin model is the first finite system which 
serves as a representation of this algebra. It can be studied in a hands-on manner, for 
instance on a computer, and one can construct the ladder operators explicitly, just 
as for SU(2) spins. 

At this point we comment on the third issue: the relation to experiments. The 
utility of the ISE-model and its Yangian doesn't lie so much in the fact that it clarifies 
experiments, waiting for an explanation, but rather, it presents a new paradigm, a 
language. We mean this in the same sense that free boson and fermion models are 
paradigms that do not occur in nature, but a class of more realistic models only differ 
from them by the addition of irrelevant perturbations; these models reside in the 
same universality class, and have essentially the the same long distance physics as 
the free fixed points. Similarly, the nearest neighbor Heisenberg model is in the same 
universality class as the ISE model, but its semions are weakly interacting, as we will 
see in chapter |. 

The remainder of this thesis is organized as follows. The first part, chapters 
[l] and H, will be concerned with the "exactly solvable" qualifier of the ISE-model, 
illustrating the richness of the underlying integrability structure. Chapters |3] and 
f| are more applied and stress the paradigmatic aspects of the model. Chapter [1] 
will review the relevant facts about the solutions of the ISE spin-model, and show 
how its Yangian symmetry emerges. In other words, we will try to construct the 
transfer matrix of this model. Chapter |2] is a bit more technical in nature. In general, 
all exactly solvable models will possess a set of commuting operators of which the 
Hamiltonian is one: the constants of the motion (in classical physics, the conserved 
quantities). The usual derivation of these quantities does not work for the ISE model, 
for the precise reason that the Yangian is an exact symmetry of the ISE Hamiltonian. 
By deforming the model through allowing the particles to leave their lattice sites, via 
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a kinetic energy term, and taking the limit where the spin potential energy dominates 
the kinetic energy we can recover these constants of the motion. 

Chapter |3] describes the dynamical correlation functions of the ISE-model in a non- 
zero magnetic field, where some knowledge of the underlying Yangian structure — the 
division of the Hilbert space into Yangian multiplets — is actually required. It builds 
on recent work based on various methods, such as supersymmetric matrix functional 
integrals [0, 0] and Jack polynomial techniques [0] from mathematics, which obtained 
these structure functions (or Greens functions in the corresponding quantum field 
theory), exactly, at at zero external field. We investigate the case h ^ 0, in which 
case the /i-field tunes the number of fractional quasiparticles in the groundstate. 

Finally, chapter |] show how a new order parameter for antiferromagnetic spin 
models [0] can be analyzed in the free spinon ISE-model. The order parameter involves 
inserting two sites into a spin chain and putting two spins in a singlet configuration 
on them. Whereas only numerical and circumstantial evidence exists for the Long 
Range Order (LRO) of this order parameter in the NNE-model, it can be approached 
analytically in the ISE-model, once we interpret the order parameter as a spinon 
propagator. Some analogies to the v = | bosonic FQHE are drawn, and we suggest 
some implications for stabilizing superconductivity in layered cuprates via singlet pair 
hopping. 
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Chapter 1 



The Inverse Square Exchange 
Spin Model and its Yangian. 

1.1 Introduction 

Solvable models play an interesting role in many-body physics and not just because of 
the mathematical satisfaction of being able to obtain the exact solution to a particular 
problem. In a world of uncontrolled approximations, they can be the touchstone 
and foundation of theories that describe "real" systems more realistically than the 
mathematical toy models themselves. The 1-D Hubbard model of spin-1/2 fermions 
on a circle illustrates this eminently [0]. 

Although nowadays scores of models have been solved exactly, many of them still 
await practical application. This is not surprising considering the mathematical tour 
de force that is often required. A consequence is that solving a model (i.e., knowing 
the eigenf unctions of that particular Hamiltonian) does not mean that we can actually 
compute physical properties if those wavefunctions are of a complicated nature. Such 
is the case for the Nearest Neighbor Exchange S = 1/2 Heisenberg-chain, for example. 
Its Hamiltonian is given by [0]: 

N 

-f^NNE = J X! • Sj+i, (1.1) 
i=l 

where Sjy + i — S±. 

A system that doesn't suffer from this problem is the model of spin-| particles on 
a one-dimensional ring, interacting antiferromagnetically through an Inverse Square 



8 



Chapter 1: The Inverse Square Exchange Spin Model and its Yangian. 



9 



Exchange interaction (ISE-model for short). It was introduced in 1988 by Haldane 
and Shastry [0, 0]. The Hamiltonian is as follows: 



ISE 



J N 

7 E d(i 



N 



2 (Sj • Sj - 



:i.2i 




where permutes the spins on sites i,j 
and d(i — j) represents the chord distance be- 
tween two points i,j on a circle that has 
been subdivided into equidistant sites: d(n) = 
sin 2 (^)| 2 . Periodic boundary conditions are 
assumed. 

Despite the long range of the force between 
the spins, this model is solvable to a much larger 
extent than its short-ranged cousins. The key 
feature is the non-transcendental structure of 
the eigenf unctions of ( |1.2| ), and the simplicity 
of the spectrum. It may seem that the particu- 
lar exchange integral in eq. (|1.2| ) makes -ffisE less suited as a paradigm for a family 
of physically applicable models (but see [0] for exceptions). However, the fact that 
the model's solution is simple, in combinations with the fractional statistics nature 
of the elementary excitations in its spectrum, make it an appropriate candidate for 
an ideal gas of anyons. 

We will start with a short introduction to the ISE-model and its integrability 
structure; precise details can be found in [0, 0]. Section IOI will explain how Hise 
was solved originally. Section L3 will analyze the symmetries of the Hamiltonian 
and provide an explanation for the exact degeneracies of its spectrum. The meth- 
ods on which section |1.3| are based, are familiar from quantum inverse scattering 
and the Yang-Baxter equation literature — thereby elucidating the relation to more 
conventional solvable models. 
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1.2 The ISE-Hamiltonian and its Spectrum 

To solve the Schrodinger equation for i^isE let us first transcribe it into a language of 
hard-core bosons on a lattice, where down-spins are considered particles and up spins 
represent the absence of these particles. The hard-core constraint refers to projecting 
out the states in the Hilbert space with more than one particle per site. With M 
particles present, S z = \N — M. The usual separation of the Sj • Sj term into a 
kinetic energy part (S^~Sj + SjS^) and a potential energy part S-S*, reduces eq. 



L2j) to [0]: 



N v s N 

Yj e o(k)pk + — d^-j) 2n i n j, (1-3) 

fc=l 71 i^j=l 

where the discrete Fourier transform of the exchange is eo(k) = ^k(k — 2ir), k = 
jfm, m — . . . N — 1. Sj = c], nj = CjCj = \ — Sj, and pt = & k C)~ is the usual 
number operator for the extended orbital with momentum k. From now on we will 
measure energy in units of (^) 2 ^. Notice that the kinetic energy is strictly quadratic 
(albeit periodically repeated in successive Brillouin zones), despite the fact that the 
particles hop on a lattice. This peculiar connection with the kinetic energy in the 
continuum is a special property of the 1/ sin 2 potential that people have attempted 
to employ before in the context of fermion lattice gauge theory [0] . This phenomenon 
raises the suspicion that the relation between the lattice model and some continuum 
model of particles on a ring may be exploited. Haldane and Shastry [0, 0, 0] found 
indeed that the set of eigenfunctions of the continuum Hamiltonian of the so-called 
Sutherland model [0]: 



M M 



'--^g gy'&W (L4) 



of M particles on a circle of circumference N are also eigenfunctions of //ise v i & 
the particle/down-spin analogy. Let us repeat this derivation since it involves a few 
transformations that we will run into frequently. 

Sutherland found the (unnormalized) groundstate for his continuum model to be 
of the simple Jastrow form: 



i<j 



sin( -7r) 

V N 1 



(1.5) 
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This is the case of bosons; for fermions we adopt Q1.5Q for xi < x 2 ■ ■ ■ < % and 
multiply with the appropriate sign for other orderings of the particles. 

As a notational aside, in the spin model on the circle — with a lattice spacing 
that we have chosen to be 1 — let us parametrize the sites with the complex numbers 
\zj = e 2 ^- 7 j . The locations of the particles (down-spins) {n%}iLi in a given spin- 
configuration denoted by the ket |ni,n 2 , . . . , n M ), are given by z nj , j — 1, . . . , M and 
< rij < N. As an example, in this complex notation the "interaction potential" 
between down-spins becomes: (% ) 2 ■ -, — 2 "' z "> . Haldane and Shastry showed that 
( p..5[ ) is an eigenstate of spin Hamiltonian ( |1.2p as well if we restrict the Xi to the 
lattice, i.e. to be integers {rii}, and pick A = 2. In other words, the spin groundstate 
is a member of the family of wavefunctions (in complex notation): 

i<j i=l i<j i 

(1.6) 

In the spin model, the absolute groundstate as function of M occurs when J = 1, M = 
N/2, and S z = 0, (at least for N even, which is what we will assume in the rest of this 
chapter). In the sector M^|we have to choose J = 1 + ^ — M to find the state 
with lowest energy for that value of S z . For what follows we choose J — 1. Here, 
and in the following we adhere to the convention that uppercase greek letters refer to 
wavefunction of the {n{\ and lowercase greek letters to functions of the {^ ni }- That 
is: is related to i/j as follows: 

^ (^l, . . . ,71m) = i>o{z ni , ■ ■ -,Z nM ). (1.7) 



As we shall see in section |1.2.2j , the excited states in the continuum are ob- 
tained by multiplying the groundstate wavefunction by a symmetric set of plane- 
wavefunctions — \1> = $? (xi . . . x M )$( 

X\ . . . Xy[ ), with momenta that satisfy simple 
Bethe-Ansatz type equations. The rationale is that the strong correlations are de- 
scribed well by the "groundstate" wavefunction and small fluctuations are added by 
the multiplying factor. In the spin case we can try the same: we write a spin wave- 
function as the product of the wavefunction (|1.6j ) and some symmetric polynomial 
in the z ni : 1 ^(ni, . . . ,n M ) = ^/ (rii, . . . ,n M )$(ni, . . . ,n M ). The crux the map- 
ping of Hise onto the continuum Sutherland Hamiltonian is the following: the spin 



1 {z n ) m is a plane wave of momentum jtj-m. 
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Schrodinger equation for ty § which is really a difference equation for $, can be 
rewritten as a differential equation in $, if the latter is of a polynomial (i.e. plane 
wave) form. This differential equation is identical to the continuum Schrodinger equa- 
tion that ...%) = |cont[^"'^j satisfies (the superscript "cont" indicates that 
the Xi are considered continuum variables). This equation has been solved in [0]. Let 
us first illustrate the difference-to-differential transformation. 



1.2.1 From the Lattice to the Continuum 

In complex notation: 



#ISE - H x + H 2 - E 7- * J \ 2 S i S j + E f v ^j\ 2 U S t S j + S i S P 



ijLj \ Z i Z j) ijLj \ Z i Z j, 



We have dropped the additive constant part in i^isE fc> r now, and we will put it back 
in at the end, by requiring that the fully polarized state | ] ■ ■ ■ ]) has zero energy. 
Now (ni, . . . , riM\Hi\§^ ) = . . . , %)(ni, . . . , nM\Hi\ty ) since Hi is diagonal 

in the particle number occupation basis. Furthermore 



(n h ...,n M \ H 2 |W ) = ^^ 



M N 



, Zi - Z: 



\2 



j=l i^rij \ * n j f 

* Q (ni, . . . , n 3 -_i, i, n j+1 , . . .)$(ni, . . . , i, n j+1 , ...). (1.9) 

We can write this as a "continuum" expression with the help of the fact that the 
Fourier transform of 1/ sin 2 is quadratic: 

E ( / n _ Z :, 2 P(Zn) = {-IMJ 2 + %Z m d Zm - ^ - 1)1 P{ Zm ) (1.10) 

if P{z) is a polynomial of degree < N in every variable. Since </>o is of has powers of 
z ni in the range 0, . . . , 2(M — 1) + 1, this means that we can only apply eq. ( |1.10|) if 
has a degree in the range 0, . . . , N — 2M in all its entries. This somewhat arbitrary 
constraint from the point of view of the spin model will become clearer later on. Then, 
with the Leibnitz rule for differentiation: the H 2 matrix elements ( |1.9|) become, 



M 



1 N 



-*(m, . . .,n M ) ( E [-^ n2 - + - i(*A) s 
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M 



N 



3=1 V 2 7 

= $(ni, . . . ,n u ){n u . . . ,n M \H 2 \^ ) + # (ni, • • • 



.n A/ |i^|$> 

(1.11] 



where H' 2 = 



s z n A. - deg($)(f - M). u; 



Zi~\-Zj 



— zcot(^7r). Remember that is a symmetric polynomial and therefore homo- 
geneous of degree deg(0) and satisfies Euler's theorem. Now using the Schrodinger 
equation: H(^ <£>) = S-($* ) and H^ = E Q (M)^ , we have H' 2 <£> = (E-E (M))$. 
This reduced Schrodinger equation for $ is identical to the equation obtained when 
\l/ ■ $ would have been considered a continuum expression and be plugged into the 
continuum limit Schrodinger equation ( |1.4| ), provided A = 2 [0]. For good measure: 
E (M) = \M(AM 2 - 1) - \M 2 N. The case A ^ 2 would correspond to an XXZ 
generalization of the isotropic spin model. We follow the solution of Sutherland for 
diagonalizing H' 2 in the next section. 



1.2.2 Eigenstates of the Sutherland Hamiltonian 

The (non-hermitean) Schrodinger equation for $ can be solved by choosing a basis 
in the space of the symmetric polynomials, such that H' 2 becomes triangular. The 
eigenvalues can then be read off from the diagonal. The obvious basis for symmetric 
polynomials are the functions 0' m ' = Y^ P£Sm IT^i Znf^ , where < m x < . . . < 
nT'M < N — 2M. Notice: deg(0' m ') = Yli rrii. The only degree of freedom left to make 
H' 2 triangular is to order the different partitions \m\. The off-diagonal action of H 2 
on a basis state 0' m ' is to produce a linear combination of other basis states that are 
labeled by partitions obtained from repeatedly "squeezing" partition \m\ [0]. We say 
that \n\ is produced by squeezing \m\ if = rik,k ^ i,j and if for some i < j: 
rii = rrii + 1, rij = rrij — 1. This is illustrated in the next figure: 



mi 


m 2 


m 3 


m 4 


m 5 


n x 
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n 3 




n s 
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Thus, if we order the states according to whether they can or cannot obtained from 
each other by repeatedly squeezing, H' 2 is by definition diagonal. This is possible 
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since this partial ordering is transitive. The eigenfunctions of H' 2 are then labeled 
by a unique partition \m\ and consist of a linear combination of the corresponding 
basis state plus all its squeezed descendants. They are known in the mathematics 
literature as Jack polynomials [0, 0, 0]. The energy of such a state is then [0]: 

N M 
E (M) + 1 1 {^2(n 2 3 -2(M + l~2j)n 3 }~(--M)J2n J = 1 ^ m j( m 3-N), (1.12) 

j 3 i=i 

where the rrij = rij + 2j — 1 are a more convenient set of variables in the range 
< rrij < N. These variables, which are distinct and cannot be consecutive either, 
are called pseudomomenta. The rrij satisfy an equation that makes contact with 
Asymptotic Bethe Ansatz methods: 

e *™ mi -N _ TT ^(rm-rrij) (1-13) 
i<j 

and 0(mj — rrij) = 7rsgn(mj — rrij). Q(m) is the quantum mechanical phase shift for 
a scattering process of two particles, with relative momentum ^Ep, that interact via 
a A(A — 1)/ sin 2 (x) potential, at A = 2. 



1.2.3 Missing States and their Origin 

The states we have found thus far are all at the top of their spin multiplets: S = S z , 
which can be seen easily from the fact that rh is a polynomial which vanishes at Zi = 
(see Appendix A). By acting with S~ we can recover the rest. This gives us a total of 
(for N even) N M ^j spin- multiplets with S = N ~£ M (we have to choose M n^s in the 
range 0, . . . , N — 2M). We know that there ought to be (^j — > (^n") [0], 

and the missing multiplets must be of non-polynomial form, that is, they correspond 
to continuum eigenfunctions with support beyond the first Brillouin zone [0], which 
gets mapped back to [— n, it) because of Umklapp processes in the lattice case. Thus 
the trick (|1.10|) ceases to be valid. 

However, by diagonalizing small systems numerically, Haldane [0] finds that all 
these "missing" states have energies already contained in the set of energies of poly- 
nomial-type states. I.e., the energy levels have enormous degeneracies beyond the 
regular SU(2) symmetry of the model (see table [TT] later on). This hints at the 
fact that there exists some hidden symmetry generated by operators commuting with 
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the Hamiltonian (but not amongst each other). Since the introduction of the model 
various operators that commute with i^isE were suggested by Inozemtsev [0] , but not 
with the intention of reproducing the degeneracies — all these operators commuted 
amongst themselves. Such an algebra should not commute with the SU(2) Casimir 
S 2 . The two Inozemtsev invariants of the third order can be linearly combined into 
an operator that commutes with S 2 : 

^^'J^LSrS.xS, (1.14) 

ijk z tj z jk z ki 

where = z,- L — Zj and the sum YJ ' excludes terms where any of the dummy indices 
coincide. We will revisit this operator in chapter [| The other linear combination: 

I' 3 = ^'{wij + w jk + w ki )Si ■ Sj x S fe 

ijk 

oc S- ^ WijS, x Sjj = S- A (1.15) 

also commutes with S and with -£/ise, but A only commutes with ifisE and doesn't 
commute with S 2 ! A is called the rapidity operator, since A 2 measures the total 
rapidity of one of the polynomial type eigenstates, i.e., Ylj m j ( n °t to be confused with 
the translation/momentum operator which determines Yij e^ mi or Ylj m j modiV). 
Now in Bethe's traditional model with nearest neighbor exchange (NNE) interaction 
this operator has an analogue: A NNE = Y]i^j9{i — j)Sj x Sj, where 6{n) is the 
Heavyside step function = lim K _»oo coth(/«(z — j)). So the NNE rapidity operator is 
seen to be a particular limit of the analytic (hyperbolic) continuation of the one in 
the trigonometric model with ISE, i.e., Zj = e~ J — > e KJ . 

The algebra generated by Sxot and Anne has been studied in the mathematics 
literature and is known as a representation of the Yangian algebra Y(sl2) (s/2 is iso- 
morphic to SU2). The Yangian, a Hopf algebra, was introduced in 1985 by Drinfel'd 
[0] in the context of the Yang-Baxter equation (YBE), and its representation the- 
ory has been worked out subsequently by Tarasov, Kulish et al. [0] and Chari and 
Pressley [0] . The set of Yangian operators does not commute with the Hamiltonian 
in the NNE model, except when iV — > 00. However in the long-range spin model 
S and Aise do commute with i^isE for any finite N. The entire representation of 
the Yangian that they generate in this case must also commute with H, i.e., for the 
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ISE case Y(sl 2 ) is a genuine symmetry, giving rise to additional degeneracies. In 
other words: this particular spin model is a reducible representation of the Yangian! 
From the Yangian representation theory we can also understand the special role of 
the polynomial wavefunctions: they are the Yangian Highest Weight States (YHWS) 
of different irreducible Yangian representations. Every one of them generates an ir- 
reducible representation space by letting S and A act on them repeatedly. All the 
resulting states will have the same energy as their parent state, but won't in general 
be polynomials. 

Every irreducible SU(2) representation is parametrized by one integer, since there 
is only one element in the Cartan algebra — for SU(N) there are N — 1 integers. For 
Y{sl2) we characterize each irrep with a polynomial, the so-called Drinfel'd polynomial 
(a polynomial is determined by the prefactor of the highest power and the positions 
of its roots; here we only care about the relative positions of the roots) 2 . In the spin- 
chain case the degree of this polynomial is N — 2M for a YHWS with M down-spins. 
The locations of its roots, which are spaced by integer multiples of some unit called 
h, are determined by the pseudomomenta {rrii}^ and vice versa. 

One of the things the representation theory of the Yangian algebra can teach us 
is the degeneracy of the multiplets: i.e., exactly which spin-multiplets occur in a 
given irrep, belonging to a particular energy level. Haldane [0] has identified a rule 
that gives the spin-content of a particular energy level, but in the next section, we 
will go through the 'Yangian derivation" [0, 0]. We will do this by constructing 
another irreducible representation in a Hilbert space of which we know the SU(2) 
spin content in section |1.5.1| . This representation is equivalent to one based on one 
of the polynomial states, as they are made to have the same Drinfel'd polynomial. 

In order to understand the representation theory we should first discuss the ab- 
stract algebraic structure underlying the Yangian, the transfer matrix and the Yang- 
Baxter equation, in |1.4| . In the case of the NNE-model the NNE representation of 
the transfer matrix (which contains Sxot, Anne and all other elements of the Yan- 
gian as coefficients in a power series) emerges naturally if one calculates its partition 
function Tr(e~^ ). We construct it specifically in |1.5| , as well as for the ISE-model. 
The transfer matrix for the latter, unfortunately doesn't have such a nice physical 
interpretation, but its representation is reducible, explaining the multiplet structure 

2 In a SU(N) spin chain a Y(sIn) representation is characterized by N — 1 polynomials, see 
chapter g. 
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Of -HlSE- 

1.3 The Transfer Matrix, its Yangian, and the 
Algebraic Bethe Ansatz 

Historically the Yangian emerged from a study of the algebra of the transfer matrix. 
Almost all exactly solvable ID models owe their integrability to the existence of a 
transfer matrix which encodes the consistency rules on scattering processes that lead 
to the Bethe Ansatz equations [0]. The transfer matrix T(u) is an operator- valued 
matrix function of a spectral parameter u. u e C and T is a 2 x 2 matrix in our 
SU(2) case (for SU(N) : N x N). Its nth Taylor-expansion coefficients in ^ contains 
essentially the operators generated by commuting n S's and A's with each other. 
Specific to SU(2) the transfer matrix is defined by 

T a \u) = 

oo i 3 

= ^ + E-^lE<^ «, 6 =1,2 (1.16) 

{01,02,03}, the three Pauli matrices, and cr = 1, satisfy Tr(<7 Q <7g) = 25 a p. These 
Pauli matrices do not act on the local spins, but are just a convenient way to label 
different linear combinations of the four components of T. The 3^ are abstract 
operators, but in the spin-chain representation they are combinations of the local 
Sj. To follow the upcoming discussion it may be helpful to take the following into 
account: 

Jo — STot 

JiocA (1.17) 

The Yangian defining commutation relations — the algebra of Sxot and A and the 
operators generated by commuting these with each other — are encoded in the Fun- 
damental Commutation Relations (FCR) for T: 



R(u - u)(l <g> T(u))(T(u) ® 1) = (T(v) <g> 1)(1 <g> T(u))R(u - v) (1.18) 
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where R is a 4x4 matrix acting on the a, b indices of both T's but not in the space 
in which their matrix elements act, ie R commutes with the J%. Here R is Yang's 
"rational" solution (whence the name Yangian) of the Yang-Baxter equation: 



R = ut + hP, 



12- 



;i.i9) 



P12 is the operator that permutes vectors in the tensor product C 2 ® C 2 : P\2{x®y) = 
y®x. h is some complex number, called the quantum parameter, which is fixed for a 
particular model. For the ISE-model, the Drinfel'd polynomial will turn out to have 
real roots, and we can normalize h to 1. For the NNE-model, complex roots occur in 
a pattern symmetric around the imaginary axis, and it is customary to choose h = i. 
From now on we will use the summation convention with a sum over Greek indices 
running from ... 3 and Latin indices running from 1 ... 3. The FCR ( |1.18| ) can be 
written out more transparently in components: 



(u - v)[T ab (u),T cd (v)} = h{T cb {v)T ad {u) - T cb {u)T ad {v)) 



-icd 



-*ad( 



;i.2o) 



This holds in general for a chain based on for sl n spins as well. Notice that for h — > 
the transfer matrix commutes with itself; this is referred to as the classical limit. For 
the particular case of a spin-- model we rewrite ( |1.20| ) in basis of a- matrices: 



{u-v)[T a {u),T\v)\ = lie ljk Ti(v),T k (u)}, (a) 
(u-v){T%u),Ti(v)] = ^e ljk (T k (v)T°(u) -T k (u)T°(v)) (b) 

= lte ljk ({T°(u)T k (v) -T°(v)T k (u)) 



;i.2i) 



{, } denotes the anticommutator. Expanded in modes this becomes: 



TO Ji 

"m+1 ' n 

P P 
"m+1) "n 



TO Ji 
u mi u n+l 

Ji jj 



^ e ijk{J n J^-, 
( 7° J k 



^^-ijk 



TO Jk' 



(b) 



[1.22) 



This is to be supplied with the "boundary condition" J® 1 = 1; J_ x = 0. In deriving 
( CT , [Q2D we have used the fact that [T^u), T u (v)} = [T»(v),T v {u)\, which fol- 
lows from adding the FCR and the one with u and v interchanged. An important 
consequence of this last relation with // = v is that the components of the transfer 
matrix commute with themselves for different values of the spectral parameter. In 
particular, if \i — v — we have [Tt(T(m)), Tr(T(t> ))] = 0. If we expand this result in 
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powers of u,v we find [J^, J°] = V mjTl , i.e., we have a set of commuting operators. 
For the NNE model, the Hamiltonian is a member of this family, and thus we obtain 
a set of constants of motion. However [Tr(T),T ab ] = only for iV — > oo, so the 
Yangian doesn't provide any extra degeneracies. Unfortunately, for the ISE model 
ifiSE cannot be a member of that family since //ise commutes with T. We will return 
to this issue in chapter Q 

Let us look at eq. (|1.22| ) for a few particular cases. When we choose m = —1 
in (|1.22| a ) we find: [J®, J n ] = so J ° commutes with the algebra and is scalar. 
Setting m = — 1 in ( |1.22f ): [Jg, Jj£] = hie^kJ^, i.e., the J are SU(2) vector operators. 
Similarly with m = in ( 1.221 °) we find that the J° are scalars. With m = in ( 1.22j 6 ) 



we have: 

i£iik[Jl Jn] = -24+i - 2/^(^0° + 1)4 - 2h4(J° n ). (1.23) 

This last equation tells us that we can generate all the vectors J n ,n > 1 once we 
specify J and J 1( plus the trace operators J^, m < n. We will get back to the latter 
shortly. 

In order for the ensuing {J n } to satisfy the FCR, Jo and Ji have to obey a set 
of consistency relations, known as the deformed Serre relations. Our candidates for 
Jo and J x viz. S and A turn out to satisfy these relations (for both ISE and NNE 
model) [0], but since we will explicitly construct a transfer matrix that satisfies ( |1.18 ) 



and has S and A as its first moments we do not have worry about these constraints. 

To fix the set we require some additional information beyond the FCR. We 
need to specify the so-called Quantum Determinant: 

Det q = T n (u)T 22 (u -h)- T 12 (u)T 21 (u - h) (1.24) 

It can be shown with the FCR that [T ab (u), Det g T(u)] = (see [0]). The quantum 
determinant commutes with the Yangian algebra. In the particular representations of 
Y(sl2) in the NNE and ISE spin chains, it is actually a c- number function times the 
identity operator (but in general it doesn't need to be, see chapter ^). From eq. ( |1.24|) 
we see that Det g (<f>(u)T(u)) = <p(u)(p(u — h)Det q T(u) so by multiplying the transfer 
matrix by a scalar function we can choose the normalization of the quantum determi- 
nant. After an expansion in powers of - this normalization fixes the constant (non- 
operator) part of the J°. We should also note that Det,(TiT 2 ) = Detg(T 1 )Det q (T 2 ). 
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1.4 Representations of the Yangian 

Like in the case of Lie algebras, finite dimensional irreducible representations of the 
Yangian are built on a Highest Weight State which is annihilated by a set of raising 
operators: { J+ = J* + iJ%}. Furthermore the YHWS are eigenvectors of J z and 
J®. The action of the lowering operators {J~ = J® — iJ%} generates the other states 
that span the representation space. This information can be encoded in the following 
equation: 

T( M )|YHWS) = f ° ] |YHWS). (1.25) 

t n (u) and t 22 (u) are complex functions of u that label the representation. Just like 
S z acting on a maximal S z HWS determines the spin of an SU(2) irrep, the diagonal 
elements of T, t 11 and t 22 , determine the particular Y^s^) irrep that we are dealing 
with. Actually, only their ratio does, since the quantum determinant constraints their 
product, //the representation is finite dimensional and irreducible, this ratio can be 
written as follows: \sz0, — ^pn?r ? where P is a polynomial. P{u) is the Drinfel'd 
polynomial. 

The other states in the Hilbert space are generated by the repeated action of 
T-{u) =T 2 \u): 

|Ai, . . . , A r > = [J T-(A,,)|YHWS) (1.26) 

i=l 

The {Aj} aren't arbitrary numbers: we want to choose them such that the states 
that they generate are an orthonormal basis for the representation. It turns out that 
choosing them to be eigenvectors of the trace of the transfer matrix T°(u) does just 
that, since T° is hermitean 3 and consequently, its eigenvectors are orthogonal (just 
like {S~) n \S z = S Tot ) must be an S z eigenstate) 4 . With the FCR (|L2l| ) this translates 
to the following set of Bethe Ansatz equations for the {\j} [0]: 

^ n (An) _ -q A n — Aj + h _ ^ n _ ^ r (1 27) 

£ 22 (A„) j^ n) A n - \j - h 

In the spin models this procedure yields maximal S z states; the remaining states follow 
by acting with S~ = Jq . Those states in the same multiplet have the same T°(u)- 
eigenvalue since S~ commutes with the trace of the transfer matrix. This procedure 



3 At least the spin chain realization is. 

4 See chapter || for a more specific discussion. 
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of constructing a representation space of the Yangian is called the Algebraic Bethe 
Ansatz [0]. 



1.5 Specific examples: NNE and ISE chains 

An important property of the solution of the FCR eq. Ql-18| ) is that if Ti(u) and 
T 2 (w) satisfy it (and matrix elements of T x commute with those of T 2 ) then so will 
Ti <S> T 2 . The simplest solution to the FCR with R-matrix ( |1.19D is: 



T(u) = (1 + ^-)a + -a ■ S = 1 + -P 01 (1.28) 

It acts in the Hilbertspace of one spin S. The second step is only valid for S — |. 
Pqi is the operator that permutes the "real" spin with the "auxiliary" spin defined 
by the matrix indices of T i.e., Poi = \ + cr • S. This representation of T is called the 
Fundamental Representation, but is referred to in the mathematics literature as the 
"evaluation homomorphism" . The Hilbert space in which it acts is just that of a single 
spin S. One can simply verify that its quantum determinant is given by for S = ~ 
and u2 ~ {h 'T(u-h) {S+1) for general S. Furthermore: t u {u) = 1 + A + JS; t 22 {u) = 
1 + ;r — -5 and P fe + ^ = which has the unique solution: 

2u u P(u) u +!±-hS n 



= TT U + - + hS - hj) (1.29) 




We see that Drinfel'd polynomial P(u) has a string of 25* zeroes, equally spaced, and 
centered around 0. 

From it we can construct a transfer matrix for N spins by tensoring it iV times, 
for N different spins: 



Tnne(m) = 1 1 



l + ftfl + tr-Si 
u yz 




;i.3o) 



Again, the second equality only holds for Si — \. This is exactly the transfer matrix of 
the NNE model. Since it doesn't commute with i^NNE, h suggests that this Yangian 
representation is irreducible, which, as we shall see later, it is in fact. Its YHWS is 
just the completely ferromagnetic state. We can expand it in powers of - to make 
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contact with the Yangian operators. We specialize to 5/2- 

00 1 / N 

t nne («) = i + />£— i£(| + <x- Si)(#% 

n=0 a \i,j=l 

Lij = h(l - 5^)9(1 - j)Pij (1.31) 

where 9(n) is the step function, that has been introduced to keep the dummy indices 
in the expansion of (|1.30Q in powers of - ordered. Pij permutes spins on sites % and 
j. This result follows from the fact that 

-fbil ' ' ' P0im = (o CT ' S«m)-f«m«m-l-^ : im-lim-2 ' ' ' ^12»1 (1-32) 



From eq. ( |1.31| ) we find by explicit calculation: J = h J2i Sj = ^S-rot an d Ji 



— 2 — S ~2~ X^j<? x Sj, as expected. 
For the ISE model there isn't such an obvious candidate for a transfer matrix. 
However, if we compare the form of the Ji Yangian generator for the NNE model with 
what we want in the ISE case: Ji oc YU<j z --z- ^ x ( wnere we have slightly rewritten 
the expression for Ajse ( |1.15| )), we see they only differ in that 9(i — j) — ► % = —3—. 
Using this replacement as an Ansatz in eq. ( |1.30| ), we find that the resulting T with 
L matrix 

1 



T = l + EU + o-'S* 



u — L 



Ln = (\-<\j )(>,,, (1.33) 

still satisfies the FCR — for arbitrary {zj} as a matter of fact [0]. The corresponding 
product form ( [L.3QD isn't so simple to find, but was eventually traced after the intro- 



duction of the so-called Dunkl operators (see chapter ^| for a more general discussion): 

7* = h £ [9 {j 9(j - i) - 9 3l 9{i - j)) Kij 

= I £K -sgn(i-j))i^ (1.34) 

where permutes ^-labels KijZj = ZiK^-, but not the spins! These Dunkl operators 
commute with one another: [71,77] = 0. Therefore, T = 1 + -^z^Poi and thus also 



N / 

n (1+ 

i=l \ 



-P 0t (1-35) 
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still satisfy the FCR (remember: % only acts on Zj's, not on spins). We turn T into a 
spin operator by replacing K iyi+ i by Pj,i+i, once ordered to the right of an expression 
(imagine the particles having dynamical z co-ordinates, as well as spins, then the 
projection is equivalent to just acting on states that are symmetric in simultaneous 
permutations of spin and z co-ordinate). Ref. [0] shows that the projected version 
of T preserves integrability and reproduces the asymptotic series ( |1.33| ). Again it is 
trivial to compute J = h^Si and Ji = —^Y,i^j w ij^i x S 3 - + h 2 ^-^-J2i^-i- At 
this point the Z{ are still arbitrary complex number, but the requirement that Jo,i 
commute with Hjse, and thus generate its symmetry algebra, restricts them to be 
equidistant on a circle. 

Representations of this version of K (s/2) are in fact reducible: we know that i?isE 
commutes with T, and we will find that there is more than one YHWS! We will now 
elucidate the decomposition. It was mentioned earlier that all the polynomial states 
translated from the Sutherland model are YHWS in the spin model. We can check 
this via eq. ( |1.25| ): those states should be annihilated by Jq and Jf and therefore 
also by T + (u). This is verified in Appendix A. They should also be eigenvectors of 
T u (u) and T 22 (u) and we will compute their eigenvalues. Potentially there could be 
non-polynomial YHWS, but since we will find that the polynomial ones and their 
descendants exhaust the Hilbert space, we can rule this option out. 

Because of the special polynomial character of the YHWSs we can write the action 
of T n (u) and T 22 {u) on them as differential operators, just as we did in diagonalizing 
F ISE eq. (PH) . It suffices to find t n (u), as t 22 (u) will follow from Det,T(u)|YHWS) = 
t u (u)t 22 (u — h) I YHWS). In general the allowed set of rational functions t u (u), t 22 (u) 
can be obtained from just knowing the quantum determinant and the fact that ^ is 
of the restrictive form P<y p^ , with polynomial P. But to identify specific Yangian 
multiplets with particular eigenvalues of Hi SE we will consider the action of T(u) 
explicitly in Appendix B. 

There we find: 

M ( h 
r»|YHWS) = J] ( 1 + ~ h{t) _ ^ 1 I VE1W -S) (1.36) 



2 



where M = y — S z is the number of down-spins in the YHWS. The Di are the 
dynamical Dunkl operators: D p = z p d p + ^7 P , obeying the same algebra as the We 
will revisit these operators in chapter 0. We can solve the eigenvalues problem of the 
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L>i very much along the same lines as the way i?isE was (partially) diagonalized. This 
is done in Appendix C. What we find there is that the eigenfunctions of the dynamical 
Dunkl eigenfunctions are most easily described in a space spanned by polynomials of 
the form Yi^Azi — Zj) Oi Z T . The eigenvalues are organized in blocks of partitions 
\n\ — n\ > ri2 • • • > tim and are labeled within the blocks by permutations of that 
partition {P 1; . . . ,P r \Pi G Sm}, t < ML I.e., for every P there is an eigenfunction 
■0p' such that its eigenvalues with respect to the {Di} are: di = ripm —P(i)+M. Any 
symmetric function of the Di — such as T 22 (u) — will be insensitive to the particular P 
and must have an r!-fold degeneracy in that block of states. Now we have to impose 
two "physical" constraints: 1) that the eigenfunctions are symmetric — down-spins are 
hard-core bosons — , and 2) they have degree in every one of the Z{ less than N — the 
hard-core constraint 3) that they vanish when two co-ordinates coincide is taken care 
of by the antisymmetric prefactor. 1) implies that of r! degenerate eigenfunctions we 
retain only one: their antisymmetric combination (to cancel the antisymmetry of the 
Jastrow prefactor). A consequence is that two can never coincide and (because of 
the M — i factor) two di cannot be adjacent. 2) restricts the rii to be in the range 
0, . . . , N — M. Actually, since the eigenfunctions have to be annihilated by T + (u), 
which contains S + , they have to vanish if one of the Zi = (see Appendix A), so we 
have rather: 1 < rij < N — M and thus 1 < dj < N — 1. 

With these additional constraints the eigenfunctions are effectively in a smaller 
subspace spanned by Yli<j( z i ~ z j) 2 Hi z % • ^{ z ii ■ ■ ■ > z m), with $ a symmetric poly- 
nomial of order < N — 2M. We recognize this form from solving eq. fll.ll ). The 
second factor Y[(zi — Zj) comes from the fact that we only retained the antisymmetric 
combination, which contains such a factor. 

We finally find that 

e\ u ) = 




nu-M^-^)) (i-37) 



8=1 



To obtain t ll (u), whose ratio to t 22 (u) determines the Drinfel'd polynomial and thus 
the particular representation, we need the quantum determinant. Since it is a scalar, 
we can obtain it by having it act on any vector, for instance on the YHWS | | *[*... f), 
corresponding to M — 0,ip — 1. In that case t 22 {u) = 1, so Det q T(u) = 1 ■ t u (u) = 
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l + Eijljis).. = 1+ (M^KeW, where |1> is the vector of length M: (1,1,..., 1) 
and 9 the matrix with elements 6ij,j ^ i and otherwise. In the case Zj = e~\ |1) 
happens to be an eigenvector of 9 with eigenvalue ^rp. Then 

, . u + h^p- A(u + h) 
DetgTM = — ^ = -L_J (1.38) 

^ { a .,hfN±l 



A(u) = Ili=i ( u + ^(2 i s the ra tio of two strings (a set zeroes of a polynomial, 

spaced by /i) of length N. With t n (M)t 22 (M - h) = t 22 {u) = ^±^, and the 

definition of the Drinfel'd polynomial P(u) we have: 

A(tf + fe) = P(u + h)f 2 (u + h)f 2 (u) = g(u + h) 
A( M ) P{u)f 2 {u)f 2 {u-h) g(u) • 1 • ] 

Since A(tt) has an iV-string of zeroes, it is easy to see that the only way to satisfy 
eq. ( |1.39|) for a polynomial g(u) is A(u) = g{u) = f 2 {u)P{u)f 2 {u — h). Note that this 



only admits solutions for P(u) if the zeroes of f 2 are not adjacent, which we know 
from the explicit calculation of the {di} they are not. The order of the Drinfel'd 
polynomial is thus iV — 2M. A graphical example of eq. (|1.39|) is shown in fig. (|1 . 1|) . 
We have absorbed the factor of ^1 into u for now. A picture like this figure is called 



X 


= zero of f 2 (u) 








= zero of f 2 (u-h) 


X 


x XxOOXxOOOX x OXx 


X 





= zero of P( u) 


1 


2 N-1 N 



Figure 1.1: A graphical solution of eq. ( 1.39 ), with d\—l, d 2 =3, d%=7, d^=12 & 
4=15 

a motif and, as we shall see shortly, the spin multiplets present in the Yangian action 
on the corresponding YHWS can be read off quite easily. 

The rule is as follows: A sequence of Q consecutive roots in the Drinfel'd poly- 
nomial is to be identified with an SU(2) spin S — ®. In the example motif we have 
therefore: (s=i)(g.(s , =|)cg)(s=|)=(5=3)e(5=2 2 )e(s=i 2 )e(s=o); thus a7 + 2- 5 + 2- 3 + l = 24- 
dimensional representation, a 24-fold degeracy of this energy level. 

Haldane [0] had already identified this state counting rule based on numerical 
studies. He introduced a nice physical description in terms of spinous, the elementary 
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excitations, where their number is N sp = N — 2M. Once we identify the di with the 
pseudomomenta rrii later on, we will see that the absolute groundstate of the model, 
which occurs for M = ~iV, S z = 0, has spinons. The lowest excited states have 
S z = 1 and thus contain 2 spinons 5 : spinons are created in pairs! In the motif picture, 
a sequence of Q consecutive Drinfel'd zeroes is to be identified with Q spinons "in 
the same orbital". To have the total spin come out right, the spinons have to have 
spin-|, and be in a totally symmetric state. This would make them bosons, except 
that the number of "states", M + 1, varies with N sp (which allows the identification 
as semions by the generalized exclusion principle, see section |1.5.2p . 

So every degenerate energy level has a YHWS (with say M down-spins), that 
corresponds to a sequence of N — 2M integers in the range . . . N. There is a spinon 
for every integer, and Q consecutive integers refer to Q spinons in a maximal S z state. 
We will now prove this statement. 



1.5.1 The Spin Content of the Yangian Irreducible Repre- 
sentations. 

We will go back to the simple fundamental transfer matrix of the NNE : Ti(u — hpi) = 
[u — hpi + |) + ha ■ Si 6 where the spectral parameter has been shifted by an amount 
hpi. We will take the tensor product of N — 2M copies of this transfer matrix, one 
for every spinon, each copy with a different shift pi. We pick the pi as the roots of 
the Drinfel'd polynomial corresponding to a particular YHWS in the ISE-model (so 
here the Pi do not correspond to the crosses in fig. but rather the circles). All 
Si are chosen to be spin-|. 

N—2M V -L 

T(//) 6d u-h Pi + - + ha-Si (1.40) 



By having T act on the state | \ \ . . . f) (this is not a state in the spin model, but 
represents all iV — 2M spinons pointing up, the YHWS). This is obviously the — or 
at least one of the — highest weight state(s) of this representation) we see that it 
has the form of eq. ( 1.25Q — just add an auxiliary spin and read off the action of T. 



The ones that don't, derive from YHWS with a higher S z and thus higher energy; see expression 



for Eo(M) in section L2 , 

6 We have changed the normalization of the fundamental T(it) by multiplying it by (u — hpi) to 
simplify the algebra. 
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t (u) = Yli(u — hpi), t (u) = Yli(u — hpi + h), so that the Drinfel'd polynomial 

reads: = IL Le "' P ^ = ~ hpi ^ which is the same Drinfel ' d 

polynomial as in the ISE model for the YHWS characterized by the same pj. The two 

representations must be isomorphic, provided that they are irreducible [0]. We will 

assume the latter is true and concern ourselves with obtaining the SU (2) content of 

eq. flOgl ) and thus QO§). 

Now we need to show that consecutive shifts, or spinons in the same orbital, only 

occur in their symmetric, maximal spin, combination. To that end we need to show 

the following fact about the tensor product: 

T(ii) = (u — hpi + | + her ■ Si) <E> (u — hp 2 + \ + her ■ S 2 ) with arbitrary values of spins 
Si, 5*2 is reducible, such that the subrepresentation contains the highest spin S\ + Si 
multiplet, only if p 2 — pi = S± + 5*2. (For p\ — p 2 = Si + 5*2 we get only the lowest 
spin | Si — S2 1 -multiplet, which we don't care about now). A new feature of Yangian 
representations is that reducibility doesn't imply full reducibility, that is: in matrix 
form the representations might look like: 



A 








B 


I rather than 








A 








B 



The span of vectors in block A is invariant under the action of the representation, 
but it's ortho-complement might not be. If p 2 — P\ = Si + S 2 , A corresponds to just 
the "highest component" of Si S 2 , that is the Sxot = Si + S2 states. Let us show 
this by letting T act in this highest spin sector: 



T(u) 



+er 



T + cr-T 

hpi) - 



u 



u 



u 



h h 
hpi + -)(u- hp 2 + -) + h 2 Si ■ S 2 



hp 2 ) + 1)S + {\{u - hp 2 



u 



hpi))A + iSi x S 2 
(1.41) 
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where A = Si — S 2 . Using the fact that Si • S 2 = SiS 2 on a highest S^ ot state, and 
the identity: iSi x S 2 = — |[Si • S 2 , A], the vector part of T becomes: 

f =\(l + u-h Pl + u- hp 2 )S + IfaSz + u-hp 2 -(u- h Pl ) - Si ■ S 2 )A (1.42) 

(S = Sxot)- Since A is a vector operator, it acting on an maximal Srot state produces 
only Srot = Si + S 2 or Si + S 2 — 1 states. To eliminate the latter and make T 
reducible, we arrange the prefactor of A to annihilate non-highest spin states. Those 
states have Si • S 2 = SiS 2 — (Si + S 2 ), so we must choose: p 2 — pi = Si + S 2 . 

The resulting transfer matrix is again of the fundamental type, but now with 
S — Si + S 2 . We will need its shift as well. To find it we must evaluate (|1.42|) for the 
special choice of p 2 — pi- The only non-trivial step is that: 

Ps=s 1 +s 2 APs=s 1 +s 2 = ( <r — <r Ps=S!+S2i (1-43) 

where Ps=s 1 +s 2 is the projector onto the highest component. This is easily proved by 
evaluating it on states of the type (S~) n \S^ = Si, S^ = S 2 ) and using that for n = 
A z = f^g and the fact that A is a vector operator. Substituting this back into the 
expression for the combined transfer matrix we have 



h 

" Sl +s 2 = (u- hpi - hSi + -) 



[u - hpi - hS 2 ) + /i(i + cr ■ S) 



;i.44) 



2 

The prefactor just changes the normalization. The new shift is pi + S 2 . Note that 
this is not symmetric in (pi,p 2 ), a fundamental property of the Yangian, where 
A® B c^B ® A in general. For our specific case Si = |, S 2 = |, p 2 — pi = 1, 
we thus build a representation containing just S = 1 states, without any S = states 
connected to them. The new shift resides between pi and p 2 . We can add a third 
spin-| by multiplying T X T 2 form the left by T 3 , to create a S — | fundamental 
representation, provided (pi + S 2 ) — p 3 — (Si + S 2 ) + S3, i.e., p 3 = pi - 1 = p 2 - 2: 
a string of 3 points. We can repeat this Q times to create a S = ® fundamental 
transfer matrix if the Pi ■ ■ - Pq are all equally spaced, and increase from the left to 
the right in the tensor product. In other words, T acts in the space of Q spins- \, but 
leaves the subspace with S^ot = 2 invariant, and in that subspace it is isomorphic to 
a single spin S = ® fundamental transfer matrix ( |1.40| ). The shift of the total spin 



y is in the center of mass of this string of {pi}- All of this is in accord with what we 
would expect from the spin S fundamental representation, in eq. ( |1.29|) . 
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Going back to ( P-.40Q , we can thus conclude that any string of Q consecutive 
Drinfel'd zeroes indeed represents a S = y multiplet. The total spin content of 
the multiplet is just the tensor product of these "string-spins". The only remaining 
potential obstacle is that the transfer matrix (|1.40 ) can be reduced further. However, 



Yangian representation theory tells us that for it to be reducible the strings of p's 
have to be in special positions [0]. Two strings are in such a special position if their 
union is again a string (they are collinear and "touch") or if one string is a proper 
subset of the other. This is obviously not the case since all strings are separated by 
at least two zeroes of the polynomials fi{u) and fz{u — h) — the crosses in fig. fll.lp . 
Similarly, we see that the NNE Yangian representation fll.30 ) is irreducible since 



all its 2S strings have the same shifts, coincide, and are therefore not proper subsets 
of each other. 

1.5.2 The Yangian, the Hamiltonian and Completeness. 

With the knowledge that the spinon description is correct we can now show that it 
is also complete [0]. For a given M (number of down-spins), we have N or b = M + 1 
orbitals available for the (N — 2M) bosonic spinons (see fig. ( |1 •![ )), which can have 
spin up or down, resulting in y 2M ^2M states. Summing over M = 0, . . . , y 

(N even) we recover indeed 2 N states. The fact that the number of orbitals available 
to the spinons changes with M and therefore with the number of spinons N sp itself, 
gives rise to the interesting issue of one dimensional fractional statistics. Since the 
number of orbitals available to the spinons goes down by one if we add two spinons 
(AN sp = — 2A(M+1)) we see that by this token, spinons interpolate halfway between 
bosons and fermions and are therefore called semions [0]. This form of statistics, which 
is based on state-counting rather than exchanging particles (as is customary in the 
two dimensional quantum Hall effect), can be applied in any number of dimensions. 
The fact that the Yangian algebra is the raising/ lowering algebra of the semionic 
spinons, means that with it we have at our disposal a well-developed mathematical 
tool for probing fractional statistics [0]. 

To connect every Yangian representation with a degenerate eigenvalue of 1/ise we 
need to rewrite the Hamiltonian's action on YHWS in terms of the Dunkl operators 
Dp Using the trick ( 1.1UQ to convert a convolution with 1/ sin 2 into derivatives we 
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Table 1.1: The spectrum of i^isE for 6 spins. The notation for the motifs in the first 
column is identical to that in fig. |1.1| . The di in column 2 are the same as the of 
eq. (|1 .12]) , which is also used to compute the energy in the last column. 



Motif 


Pseudomomenta 


Drinfel'd Roots 


Spin Content 


Degeneracy 


E in units of 


123456 


{di} 


{Pi} 






/ 2tt \2 JJ S 

V JV / 7T 


XxXxXx 


1,3,5 







1 


-9.5 


OXxXxO 


2,4 


1,6 


i © | = 1 ©0 


4 


-8.0 


OOXxXx 


3,5 


1.2 


1 


3 


-7.0 


XxXxOO 


1,3 


5,6 


1 


3 


-7.0 


OXxOXx 


2,5 


1,4 


1 © | = 1 © o 


4 


-6.5 


XxOXxO 


1,4 


3,6 


4 


-6.5 


XxOOXx 


1,5 


3,4 


l 


3 


-5.0 


OOXxOO 


3 


1,2,5,6 


1 ® 1 = 2 © 1 © 


9 


-4.5 


OXxOOO 


2 


1,4,5,6 


\ © f = 2 © 1 


8 


-4.0 


OOOXxO 


4 


1,2,3,6 


f © | = 2 © 1 


8 


-4.0 


XxOOOO 


1 


3,4,5,6 


2 


5 


-2.5 


OOOOXx 


5 


1,2,3,4 


2 


5 


-2.5 


000000 




1,2,3,4,5,6 


3 


7 

+ 64 


0.0 



have (in units of — ): 

tfisE = IE7f^(^-i) = EiA(A-iV) 

= EH^-^) d i e {!,...,#- 1}, i = l,...,M (1.45) 

i 

which reproduces the result from the Sutherland model diagonalization if we identify 
the di with the pseudomomenta r/v Knowing the spectrum of i^isE as well as the 
degeneracies of the energy levels, one can obtain the thermodynamics [0]. In table 
[TTT] we have listed the entire spectrum and its degeneracies for iV = 6 spins. 



1.6 Conclusion 

We have found that in the ISE model we can separate the Hilbertspace of 2 N states 
into groups labeled by a sequence of M integers {di} in the range 1, . . . , N — 1 and 
M = 0, . . . , y. For every sequence or set of pseudomomenta there is one spin state 
which is of a special polynomial character, called a YHWS. The energy of such a 
state is (in units of (^r)— ) | Si di{di — N). The rest of the states in that group are 
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generated by the repeated action of the Yangian lowering operator T~ (it) on it. Those 
states are all degenerate with their YHWS. The motif picture of the pseudomomenta 
interprets the Yangian multiplet as N — 2M bosonic spin-| spinons in M + 1 orbitals. 
In the YHWS, these spinons are all fully polarized (maximal spin), and the subsequent 
action of the Yangian rotates these spins and lowers the polarization. 

For the NNE Heisenberg model, there is only one YHWS for finite iV: | TT • • • t) 
and the spinons picture is only approximately correct (see chapter |3|). 
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1.7 Appendix 

1.7.1 Appendix A: Sutherland wavefunctions are YHWS. 

In this Appendix we will show that the subset of eigenstates of i^isE that derive 
from the continuum wavefunctions of the Sutherland model are indeed YHWS. That 
is to say: they have to be annihilated by both S + = Jq and or equivalently 
Af SE . The Sutherland states are of the form ^/(rii, . . . ,n M ) = ip(z ni , . . . , z nM ) = 
Ufij{z ni - z U] ) 2 nfii z nt ■ <f){z ni ,. . . , z nM ), where <p is a symmetric polynomial, such 
that ip is of degree less than A" in every one of the z ni . Then the action of S + (which 
reduces M by one) is, using the symmetry of the wavefunction: 

(S + ^f)(n l ,...,n M ^ 1 ) = J2^( n u - ■ ■ , n M-i,j) =^^(4 1 ,--,4 M _ 1 ,e^) 

3=1 3=1 

= ^(z ni ,...,z nM _ 1 ,0) = 0. (1.46) 

The last step holds because the sum over j is just a Fourier transform at momentum 
on the last entry of ip, which picks out the (vanishing) constant part in ip. 

The second Yangian raising operator is A + oc Y^i^jWijSf . In a basis of down- 
spins S* is given by: \ - Y,f=\ 8i, nj ■ Then: 

N I N \ 

(A + ^)(ni,...,n M -i) = | E E w v ' *( n i> • • • ,n M -i,i) 

iHnp] WO J 
M N 

~E E w mk ij(z ni ,...,z nM _ 1 ,z i ) (1.47) 

The first term vanishes because the the sum over the odd function Wij is zero. We 
can simplify the second term, involving the convolution with w^, with the following 
identity, inspired by eq. ( |1.10| ): 



N 

£ WyPizt) = NP(z 3 ) - 2z j P'{z j ) - NP(0) (1.48) 

i=l 

[P 1 is the derivative of P) if P is a polynomial of degree less than N. The z = term 
vanishes as with S + ^, and the other two are zero as well since \I/ has a double zero 
when two of its arguments coincide. 
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1.7.2 Appendix B: T 22 {u) Acting on a YHWS. 

We will follow the derivation in [0]. In a basis of states {|ni, . . . , tim)} where the 
label the positions of the down-spins, we consider the action of T 22 {u) (we set h — 1 
during the calculation): 

oo i N 

T 22 (u) = 1 + £ — £ (L% Xf (1.49) 

n=0 U i,j=l 

with Lki = h(l — 5ki)0kiPku &ki = ^ , and the z's lie on equidistant points on 
a circle. X 22 is the operator that projects out the states with a down-spin on site 
j. We let T 22 {u) act on a polynomial type state if)(z ni , . . . ,z nM ) of which the only 
constraint is that it is symmetric and vanishes when one of the Z{ become identical 
or equal to one of the other Zj. Thus the powers of the z\ in every term are at least 



one. This corresponds to ip being a maximal S z state (see Appendix |1.7.1|) . Under 



these circumstances we have (eq. [L.48|) : Z)fc(^n) ^fen-P(^fc) = ZkP'{zk) ^-P(zk), if P 



is of degree < N — 1. Evaluating the n = term in (|1.49|) is trivial. The n — 1 term 
for fixed i Yfj OtjPijX 22 gives: 

^ {1 \i\n 2 , . . . ,n M ) = (i,n 2 ,. . . ,n M \Y^' OijPijX 22 \ip) 



N 



'Qij^ti, n 2 ,..., n M ) + 9in^{h n%,..., n M ) 



i=i 3 

= {Zidi - : t— ) < 0(^kn 2 ) • • -i z n M ) + '^rw*(*> n 2, ■ ■ ■ > %) 

2 J 

(1.50) 

We have used the fact that ip is symmetric in its entries. The prime on the summation 
symbols indicates that the terms with diagonal elements of 9 are absent. The second 
term comes from the exchange from two down-spins and the first is the result of an up- 
and a down-spin trading places. Notice that the wavefunction ty 1 is only symmetric 
in entries 2 through M, and has one down-spin fixed at the (parametric) site i. 

For the n = 2 term we have to convolve this set of functions with L resulting in : 

^ {2) (i\n 2 , . . . ,n M ) = (i,n 2 ,...,n M \^2'L ij L jk Xl 2 \ip) 

jk 

M 

j 1=2 
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h((D 1 - ^±1)*W ) (i\n 2 , ...,n M ). (1.51) 



Dj = zjdi + J2k(^j) Kjk- The {Kij} are the same operators that are used in the defi- 
nition of the Dunkl operators (|1.34j) . is obviously symmetric in entries 2, . . . , M 
as well and we can now do all orders inductively. Finally, we have to do the sum 
over all locations of the spin "z" in ( |1.49| ). Since the down-spins are indistinguishable 



and now unrestricted, we are allowed to symmetrize ^ n '(rii, . . . , Um) (calling final 
summation index i now ni). This is done by applying J2fLi K\j where K\\ = 1. It is 
easy to verify that K^Dj = DiK^ so with = \1> (0) : 

N in M Mi 
n=0 j=l j=l J 



It can be shown that the summation (|1.52| ) can be brought into one fraction by 



exploiting a special case of the Calogero-Sutherland model with spin (see chapter 
eq. (|2T7|)) by acting on a state in which all particles have their spin up. 



T 2 \u) = n ( 1 + _ h{ ^_^ I (1-53) 
with 

Dj = Zjdj + '°jkK jk - K 3k 

k k(<j) 

= Z A + 5 £ ( w Jk - sgn(j - k)) K jk (1.54) 





1.7.3 Appendix C: Dunkl Operators and their Eigenfunc- 
tions. 



The Dunkl operators from eqs. ( |1.36| , |1.54| ) are given for arbitrary odd A by D{ = 
Zidi + | ( w ij — sgn(i — j))Kij. We have introduced the free parameter A so 



that we can accommodate A = 1 in chapter |1 . 1| , and A — > oo in chapter ^. Even with 



arbitrary A the Di still satisfy the Hecke algebra relations (see chapter |2|). We will 
follow the derivation of [0]. 

We expect the eigenstates to be polynomials, but we have to select some subset of 
polynomials in order for the z ^ Zj term in Di to still produce a polynomial. Therefore 
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we choose the form ip(z ly z M ) = i>o(zi, z M )<f>(zi, z M ) where = Ui<j(z l 



Zj) to cancel the pole when two z's come together, and 4>(zi, . . . , zm) = Ili=i(^. 



The action of D, on 6 is then 



A(^O0) 



A 



+ 2 E sgn(i - - 1) + — ^— + 2 2— | 



= Ac 



;i.55) 



since | X^j 'sgn(? — j) = « — ^j-^. Now, for rij 7^ (if = rij the result is obviously 
0): 



sgn(i - j) 



^ni-rij-l m-k rij+k. 
rij-rii-l rii+k nj-k _ 



rii > rij 



z 



Za 



, (1.56) 



we see that the action of one of the non-diagonal terms is to produce a series of 
squeezed polynomials, in the sense that the degree remains the same, but the powers 
get closer (the high ones decrease, and the low ones increase). Inspired by the solu- 
tion of the Sutherland model we can attempt to define a hierarchy on the partitions 
M = { n i > n 2 ■ ■ ■ > u m} such that L>i becomes triangular, and we can read off 
its eigenvalues from the diagonal. In anticipation of what is to come, we group the 
polynomials into blocks which have the same set of powers {rii} but differ in their 
ordering. For two partitions \n\ and \m\ we say that \n\ > \m\ if, counting from the 
left, for the first unequal pair {n^mi) holds > mj, e.g. {7,4,3,2} > {6,5,5,5}. 
"Squeezing" obviously only produces states in lower 
partitions. The sgn(i — f){Kij — 1) term just connects 
states in the same partition. So Z)j is block-triangular 
in this basis, and we only need to make the blocks 
on the diagonal, which contain only matrix elements 
of states of the same partition, triangular. In such a 
block, states are all labeled by permutations P E Sm- 
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Using ( |1.55 ) and (|1.56|) and forgetting terms with the 



wrong partition, we have 



+z™ P(l) z™ P(j) [sgn(n P{i) - n P{j) ) - sgn(i - j)) 



+Zi P{]) z j P{1) [sgn(n P (j) - n P{j) ) + sgn(i - j)] (1.57) 



The third, non-diagonal, term is non-zero only when % > j and < n P ^ or i < j 
and Hp(j) > n P (jy, in either case, when ordered by i,j the power on the right is always 
goes down in the application of D[. So if the permutations are ordered by looking 
at the last non-identical pair nqn) > n P u-\, implying that |rig| > \n P \, D[ is also 
triangular inside the blocks. The eigenvalues on the diagonal read: 

4 n '' P = n P (i) s M n P(i) ~ n P(j)) s g n (* - j) + x ( M r 1+l - M 2 1 ) 

= n m + \(M-P(i)). (1.58) 

To evaluate the sum over n P ^ we first pick P = 1 and then apply the permutation 
P, which is trivial since Y^j ' is invariant under permutations. Notice that the set of 
eigenvalues {<i' n '' P }f£i for fixed M,P 6 Sm is the same for all Di. 

If two of the powers in the expression for are identical ( |1.58| ) is no longer valid, 



but in our case <fi will be taken to be antisymmetric so this will not cause problems. 
Note that the highest power occurring in the actual eigenstate labeled by partition 
| ?7. | is max{rij} + M — 1. 
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Chapter 2 



Integrals of motion of the ISE 

Model 



2.1 Introduction 

In the previous chapter we presented two different spin-model representations of the 
Yangian Algebra, Tnne(^) , eq. ( |1.30| ), and Tise(w), eqs. ( |1.33| ), ( |1.35| ). Via the 



algebraic Bethe Ansatz eq. (|1.26| ) we can construct a 2 Ar -dimensional basis of the 
Hilbert space. For T NNE we need only one HWS: | ft • • • |), but for T ISE , being 
reducible, we needed several YHWS, viz. the set of polynomial type states ( |1.25| ), 
( |1.36Q (which are are also, miraculously, eigenstates of H\se)- 



In general, in choosing a different basis for the Hilbert space we haven't said 
anything about the action of a particular Hamiltonian in that space. This is where 
the power of the transfer matrix method comes in: the specific way in which the ABA 
basis states were constructed was by forcing them to be of the form ( |1.26|) , with Aj's 



such that they were eigenstates of T°(u) V u . This automatically leads to Bethe Ansatz 
eqs. ( |1.27|) . This somewhat unintuitive constraint becomes imperative once we realize 



that i^NNE is "contained" in T°(u) in the following sense: as elucidated in the previous 
chapter [T°(u),T°(v)} = <J=> [J^, J°] = 0V m ,„. The (hermitean) J° form a set of 
commuting operators. One particular combination of them: {ln(T°(-u))} oc 
i^NNE- The ABA basis is thus a set of eigenvectors that simultaneously diagonalize 
i^NNE and the other (combinations of the) J°. Since all J° commute with -£/nne 
they are conserved and therefore called constants of motion or invariants. For Tnne 
total spin, S 2 , is in that set, as is is P123...V, the momentum/translation operator, as 
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expected. In short, diagonalizing the NNE Hamiltonian is taken care of by eq. ( |1.27 ) 



with the added bonus of obtaining the constants of the motion and their eigenvalues. 

For the ISE model we can still construct the ABA basis states, which are also 
orthogonal, but from the point of view of diagonalizing Hise this achieves little since 
[ifiSE, Tise(m)] = 0. Since it commutes with Tise, -^ise cannot be obtained from 
Tr(T ISE (w)) since [Tr(T(u)), T(t>)] ^ 0! i/isE must have a different origin. Also the 
question of additional constants of motion in the ISE model isn't resolved. We already 
mentioned one invariant: if 3 in eq. (|1.14j ), conjectured in [0], which, as it turns out, 



also commutes with ifisE- Based on the same work, an additional invariant was 
conjectured [0], but this was not the result of a systematic algorithm. 

Minahan and Fowler [0] and Sutherland and Shastry [0] introduced sets of invari- 
ants that commute with the Hamiltonian, based on operators introduced by Poly- 
chronakos [0]. However, the generating functions for these sets are essentially the 
trace of the transfer matrix and thus contain only elements of the Yangian algebra 1 . 

In this chapter we will construct the rest of the set {H n } of extensive operators 
that commute among themselves and with the Yangian. In order to do this we have to 
consider a more general dynamical model in which the particles are allowed to move 
along the ring: the Calogero Sutherland model (CSM) with an internal degree of 
freedom. This model, which has been studied in [0, 0] has the following Hamiltonian: 

N / Q \ 2 

M -V- {h+p,i) ^- (21) 

Pij permutes the spin of particles rather than sites i and j. When h —>■ 00 or equiv- 
alently h — > 0, the particles 'freeze' into their classical equilibrium positions, and 
barring some subtleties we recover the spin Hamiltonian H 2 . 

The reason for this diversion through the dynamical model to obtain the constants 
of motion is the following: the so-called quantum determinant of the transfer matrix, 
which commutes with the Yangian algebra and therefore a natural candidate for the 
generating function of the constants of motion, happens to be c-number in the spin 
model (i.e. when h — > 0), as we will verify again below. But in the general dynamical 
model, this is not the case, and by carefully taking the limit h —>■ 00 we can isolate a 



lr rhe authors of Ref.[0] claim to have found the Hamiltonian H2 in their third order invariant, 
but in reproducing their algebra we did not find any such term; in fact we only recovered Yangian 
operators. 
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generating function for the {H n }. 

We will first review briefly the integrability structure of the CSM which bears 



great resemblance to that of the ISE model. Then in section 2.3 we will construct the 



constants of motion. Finally section [2.4| will show us how to obtain the eigenvalues 
of all the constants of motion. Unless noted otherwise, we will use general SU(p), 
p > 2 spins (i.e. spins that can have p values rather than 2), since most results in this 
chapter generalize away form the SU(2) case quite easily. 



2.2 Integrability in the Calogero Sutherland 
Model 

Let us first review the role of the Yangian algebra in the dynamical model. A more 
extensive treatment can be found in [0, 0]. The integrability of the CSM is also 
based on the existence of a transfer matrix T ab (u) that commutes with the CSM- 
Hamiltonian: 



OO 7 

T ab (u) = 5 ab + Y— -T" b 

n=0 " 

N 
i,i=l 

Lij = SijZjd Zj + (1 - 5ij)h0ijPij, (2.2) 

where Xj b , a,b = 1, . . . ,p acts as \a)(b\ on the spin of particle i, and 6^- = - z ^ z , as 
usual. This transfer matrix satisfies the Yang-Baxter equation as well: 

Rooiu - v)T (u)T°'(v) = T°»T° (u)R 00 ,(u - v), (2.3) 

with Rw{u) = u + hPoo, and T°(u) = T(u) <g> 1 and T°'(w) = 1 ® T(u). For the 
purposes of this chapter we will deal with another form of the same transfer matrix. 
We reiterate the representation of the so-called Z?wnA;/-operators of eq. (|1.36|) : 



h 

D,i = TiZid Zi + 7i = tiZid z . + - Y ( w ij " s g n (^ ~ j)) K ij 

wij = Z ~^ 3 - (2-4) 

Zi Zj 
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We have made the scale of Di explicit by putting a factor % ^ in front of z p d p . 
The idea is that eventually h — * 0. Kij is the operator that permutes the spatial 
co-ordinates of particles i and j. These Dunkl-operators commute: 

[A,4-] = 0, (2.5) 
but are not covariant under permutations: 



[K ili+1 ,D k ] = if fc^ + l 

K^+iDi - D i+1 K iji+1 = h, (2.6) 

defining a so-called degenerate affine Hecke Algebra. In terms of these Dunkl operators 
we can define a transfer matrix that also obeys the Yang-Baxter equation: 

TV) = (l + • • ■ (l + ^f] • (2.7) 

V u-DiJ \ u-D N J 

It satisfies eq. ( |2.3| ) trivially since the {Di} commute amongst themselves and com- 
mute with the Poj, since the latter only act on spin degrees of freedom; furthermore 
V ^ u-d- ) * s ^ e e l emen ta r y transfer matrix with spectral parameter u — Di. To 
retrieve T°(u) from T°(u) we have to apply a projection II, familiar from chapter 
|I] to T°. It replaces every occurrence of K iti+1 by once ordered to the right 

of an expression (this is equivalent to having the unprojected operator act on wave- 
functions that are symmetric under simultaneous permutations of spin- and spatial 
co-ordinates) [0]. It is now possible to validate eq. ( |1.53|) in Appendix B of chapter 



[T]: it follows from (|2.7|) and (|2.2| ) by having it act on a state where all dynamical 
particles have the same spin (say up). We will drop the "0" superscript on T(u) from 
here on. 



2.3 Construction of the Constants of Motion 

It was shown in [0], that this dynamical transfer matrix commutes with the CSM 
Hamiltonian. As does the transfer matrix, the quantum determinant must commute 
with Hise- It is given in its full SU(p) glory by: 

Det,(T(w))= e(a)T lai (u)T 2a2 (u-h)---T p(7p (u-h(p-l)). (2.8) 
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It satisfies [T(-u), Det 9 (T(w))] = and has been computed in [0] as: 

Det 9 (T( U )) = n Det, t «))n = n \ ; n. (2.9) 
where (making the dependence on h explicit): 

A(u,h) = H(u-D t (h))- (2-10) 

i 

Now obviously: 

[A(u,h),f(v,h)] = 0. (2.11) 

This holds since the Z)j's commute with each other and the Pq/s . The projector 
doesn't get in the way since a product of projections is the projection of the product— 
both T(y) and A{u) are symmetric under simultaneous permutation of spin- and 
spatial co-ordinates [0]. Since we know the eigenvalues of the commuting set Di form 
Appendix C of chapter [l], the eigenvalues of A(u) are also known: for every partition 
\n\ there is an eigenvalue: 

N AT _l_ I 

*W(«) = U(u - hn 3 - h(j - —I-)) (2.12) 

i=i 

We notice that as h — > 0, i.e. in the limit of the ISE model, all eigenvalues become 
identical and A(u, 0) is a multiple of the identity operator (compare eq. ( |1.38| ). Thus 



no non-trivial constants of motion are contained in A(u, 0). Nevertheless let us study 
( CT ) for small h. Writing f(v, %) = E„ h n T n (v); A(u, %) = En h n A n (u): 



= [f(v,h),A(u,h)} 

= [f (t;),Ao( M )] + ^([fo(t;) ) A 1 ( M )] + [f 1 (t;),Ao( M )])+0(/i 2 ). (2.13) 

The 0(h ) term is trivially 0. The rest of this letter will focus on the vanishing of 
the 0(h) term. As we shall show below: [Ti(t> ), Ao(u)] = 0. Therefore we have the 
important result: [T (f ), Ai(it)] = 0, i.e. the 0{K) term in A(u,h) commutes with 
the transfer matrix (and therefore the Yangian) of the ISE spin model. Furthermore 
it will also become apparent that 
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[A 1 (u),A 1 (u')] = 0. (2.14) 

So we can take Ai(u) to be the generating function of the constants of motion of the 
ISE model! In order to establish these results we first need to prove the following 
corollary: Zid Zi Ao(u) = when evaluated with the particles at their equilibrium 
positions, i.e. Zj = exp f^pj)- 

From eq. ( p.lOj ) we have A (u) = Hi(u — ji). Since we know that A (u) is scalar 
we can evaluate it by having it act on any convenient state, e.g. the one where all 
particles have identical values for their internal degree of freedom (for p = 2 we would 
say: all spins pointing up). I.e. the permutations reduce to 1. This has been done in 
[0]: 



A (u) 



det(u - 6) 
hz 



Z{ Zj 



1-5,: 



(2.15) 



Then using d x detA(x) = Tr[A 1 d x A(x)]det[A(x)] we have: 



-d Zi A (u) = Tr 



u-Q 1 



det[w - 0] 



J2 w" (n+1) Tr [Q n d Zi Q] det[u - 0]. 

n=0 



(2.16) 



Now evaluate the trace in a basis where is diagonal. This can clearly be done for 
— > and Zj = exp (^pj)- has eigenvectors ip n , where (ip n )j = -^=exp (jjfjn), 
with eigenvalue h f^y^ — n\. Then 



A' 



/N+l \ m 

Tr[0 m ,9 2i 0] = J2Md Zi e\i>n)h m - n 

n=l ^ 1 



Using that 



(0*e) 



jk 



(zi-zj) 2 ( Z jfiki Z k 5ij) j 7^ 



k 



otherwise 



(2.17) 



(2.18) 



we find: 
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WzM^n) = E Tf—T^ —2 



p (^t) ( e w p — e at J 
ih ^ sin fe) 

^£^tr a (2 - 19) 



Therefore d z .A Q (u) = at Zj = exp 

From expanding A (it, Ti) to 0(h) in ( |2.10 ) we have 



N I i— 1 \ /AT \ 

= e nc« - 7.) *a n (« - 7,o • (2.20) 

i=l \j=l / \j=i+l ) 

Then, with the corollary and the fact that [7i,7j] = for all i,j (the Dunkl algebra 
( |2.6|) is satisfied for h = as well): 



N (i-l \ /at 

[Ai(w),Ao(w)] = E II(«-7i) M„Ao( M )] II(«-7i) 

i=i \j=i y \j=i+i 

AT /i-l \ /AT 

= E IK«"7i) {^Ao(u)} II (« - 7i) ) = (2.21) 

i=l \j=i / \j'=i+l 

With this result we find: 



1 1 . 1 * , . 1 



ft^ Ai(u)- + Oihf 

A(u,h) A Q (u) A (u) J A (u) 

1 . 1 * , . „,, 2 



AM 



h ^Ai(«) + C(r). (2.22) 

Ao(u) 



This can be checked by multiplying the LHS and RHS by A (it, K). Now let us expand 
T(v,h) to 0(7i): 



AT / hP \ 1 ^ 

T(«,ft) = n =F^n(«-A+^ 

<=i V v-DiJ A(v,h) i= i 



1 v 



x ,11(^-7, + ^) + 

A (u) i=1 
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Ao(v)J 'th 



[v-% + hP ( 



0/ J 



N i-l 



N 



-XTT E II(« " 7i + hP oj ) zA, n (« - 7i + Wfy 
t (i;)+nfi(T;) + O(n 2 ) 




It is now obvious how [T 1 (f ), A (m)] = 0. The contribution from the first term in 
curly brackets in ( j2.23| ) vanishes by virtue of eq. ( 2.21 ). As far as the second term is 
concerned: the 7, commute amongst each other and with the P Q j, and [Zid Zi , A (u)} = 
as we showed before. 

So far we have established that Ai(u) respects the Yangian symmetry, but we also 
need to show that it is a good generator of constants of motion in that it commutes with 
itself at a different value of the parameter u: [Ai(it), A 1 (w')] = 0. It will be enough 
to prove this on the space spanned by the Yangian Highest Weight states, the highest 
weight states of a representation of the Yangian algebra. All other states in the model 
are generated by acting on these YHWS with the elements of the Yangian algebra, i.e. 
the transfer matrix. Since Ai(tt) commutes with T (t>) for any u, [Ax(u), Ai(tt')] = 
will therefore also hold on the non-highest weight states. First of all we should note 
that Ai(tt) and Ai(w') leave the space of YHWS invariant. This follows from the fact 
that all such states |T) are annihilated by TQ b {v) with a < b (see ref. [0]). But since 
T£ b (v) commutes with Ai(«), T afc (Ai(u)|r)) will also be for a < b. 

The proof that Ai^and Ai(w') commute hinges on the existence of an operator 
in the SU(p) ISE model that commutes with both these Ai's and is non- degenerate 
in the space of YHWS. To find such an operator we need to review quickly how the 
action of the transfer matrix on a YHWS generalizes from 577(2), eq. ( 1.25 ) to SU(p). 
In that case [0, 0]: 



f (m) I YHWS) 



V 



t n (u) 





tP p (u) J 



I YHWS), 



(2.24) 



where * is the set of lowering operators in the lower triangle of T {T^w)}^-, and 



t aa (u), a 



,p are scalar functions. Rather than one Drinfel'd polynomial, the 
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X = zeroes of P t {u) 
YV = zeroes of P 2 (u), P 2 (u-h) 
Z,Z; Z = zeroes of P 3 (u), P 3 (u- h ), P 3 (u-2h) ^ ? ^ * ^ ( ^ J> <? <Y v) (Z * ,) (Y y) (X) 

• = zeroes of A(u) 12... ...16 

* = pseudomomenta 



Figure 2.1: Example of a string motif for SU(3), for N = 16 spins. 

p > 2 case has p — 1 polynomials Pk{u) describing the HWS, for the p — 1 elements 
of the Cartan algebra. They are defined through: 

n'^=^ i;t =l p-l, (2.25) 

P k (u) tPP(u) 



k=j 



and t pp (u) is again of the form as in eq. ( |1.37|) . The fact that the eigenvalues 

of T on a YHWS have to be of the form ( |2.25|) severely restrict where the roots of 
the polynomials can lie. E.g. given that the eigenvalue of the quantum determinant 
belonging to this YHWS is given by ^7^p ( |1.38| ), with: 

A(u) = f[ t "(u-h(j-l)) 

N 

= H(u + h(^-k)) (2.26) 

k=l 

we find that the set of integers from 1 to N (modulo the factor ^rp) — the zeroes of 
A(u) — must be subdivided into strings of length no greater than p. Every root of the 
Drinfel'd polynomial Pk{u) labels the leftmost point in a string of k zeroes of A(u). 
These zeroes come from Pk(u), Pk(u — h) , . . . and Pkiu — h(k — 1)) respectively. This 



is most easily explained in a figure (see fig. |2.1|) . 

It is clear from knowledge of the action of T on a particular YHWS, i.e. the set 
{t kk (u)} p k =1 , we can uniquely determine the appropriate set of Drinfel'd polynomials, 
since the solution of an equation of the type = P ^ + ^ , for polynomials f,g 
and P is unique (up to a constant). As a matter of fact, just knowing * p J"s = 
Fl p^tt)---p'~l(M) ffc ' > we can recons t r uct all Pi(u) since every root of f p (u) has to cause p 
consecutive zeroes in A(u). 2 So if * p J"j uniquely determines the YHWS, T 11 is the 
nondegenerate operator that we are looking for. 



2 For the sake of completeness: to relate a particular set of Drinfel'd zeroes to the pseudomomenta 
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We can now prove that [A(u),A(u')] = 0. If |T) is a YHWS with motif T then 
Ai(u)|r) and Ai(w')|r) both are scalar multiples of |T) since they have the same 
T n (u)-eigenvalue ([T n (v), A^u)} = [T^(v), A^u')} = 0). So in this YHWS-space 
both Ai(u) and Ai(vf) are diagonal and thus commute. 

In the remaining part we will reproduce the integrals of motion that have al- 
ready been found [0, 0], and point out some subtleties in their construction. As 
is customary for the Heisenberg model with nearest neighbor exchange we take 
Ti(u) = ^-ln(Ax(u)) rather than Ai(u) to be the generating function for the in- 
tegrals of motion. In the NNE model this need arises, in order for -£Tnne to be the 
local, i.e. to just couple nearest neighbors. An added advantage is that the spectrum 
of these operators will be additive (in both NNE and ISE models), reaffirming that the 
parameters that characterize the eigenstates — be it the pseudomomenta in the ISE 
model or the rapidities in the NNE case — describe single quasiparticle excitations. 

When expanded in powers of it, Ti(u) reads: 

N 1 

r x («) = nE — — - — -utoO(h) 

fr[u- hzid Zi - 7i 

OO ( N 71-1 

= E u-^ nE E^A^r^n 

n=0 [ i=l p=0 

oo 



Y,u- {n+1) H n , (2.27) 



n=0 



where we have reinserted the projection operator that turns Kij —>■ Pij, when ordered 
to the right of an expression. We have worked out the first few H's. With = Zi — Zj 
we have: 



N 



Hi 



i 

E 



i=i 



E'^T^^ + ^E'a-K) 2 )^ 

ijk z ij z jk z ki ij 



that go into the expression for energies eq. ( 2.28 ) and ( 1.45 ): they are the roots of the polynomial 



AM 



[[ v j=1 p j^-h(j-i)) 



This is essentially A(u) with the last root in every string left out. 
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ijk Z V Z jk z kt 

H 4 = J2' (PiJM-1)- _ 

ijkl Z ij Z jh z kl z h ij 

A prime on the summation symbol indicates that the sum should be restricted to 
distinct summation-indices. To compute the previous expressions, we normal ordered 
the Zid Zi to the right in eq. (|2.27|) and thenput Zj = exp {^jfj\ The identity WijWjk + 
WjkWki + WkiWij = — 1 which lies at the heart of the integrability of these ^-models 
is very useful in the reduction of these expressions. P indicates the total momentum, 
or degree of the polynomial YHWS, and we will discuss its interpretation later on. 
Notice that fortunately i^isE is a member of this set, as i?2- Notice the absence of 
Yangian operators, as well as terms containing both permutations and derivatives. 
Although one would expect this, based on the fact that, physically, the dynamical 
degrees of freedom freeze out when h — > 0, it still seems somewhat accidental that 
they cancel in an explicit computation. The expressions in eq. ( [2.28| ) can be seen 
to coincide with those reported previously [0], 3 lending credibility to this way of 
deriving the integrals of motion. Unfortunately for large n it becomes prohibitively 
complicated to compute H n . 



2.4 The Spectrum of the Invariants. 

We have an alternative way to verify the validity of these constants of motion as well. 
We will proceed to compute the eigenvalues of the operators H n and compare these to 
the 'rapidity' description of the eigenvalues in [0]. We will constrain ourselves to the 
SU(2) case to simplify the algebra. The method that we have in mind is to produce 
the eigenvalues of A(w, Ti) of the dynamical model and take the appropriate limit 
h — > 0. This method, called the Asymptotic Bethe Ansatz is designed to work in the 
low density limit. In that case, between collisions, the particles are in free plane wave 
states. Then, by taking one particle and bringing it all the way around the circle and 
scattering it through the other N — 1 particles, we can constrain their allowed "free" 
momenta via the phase shifts picked up in every collision process. The spectrum of the 

3 We should point out a correction in eq. (7) of Ref. [0] where should be replaced by +^H2- 

This changes the invariant in a harmless manner, but this is relevant for comparing the eigenvalues 
of the operators in that article and the ones that we will find later on. 
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model then follows from adding the free kinetic energies. For most integrable models 
it has been shown that, remarkably, this procedure is actually exact, or at least can 
be made exact with small modifications [0]. From eq. ( |1.13| ) we know that the ISE 
model is in this category, as are the NNE model and the continuum Sutherland model 
( |1.4|) . For the CSM, integrability has been confirmed [0], so one might expect that 
in that case the asymptotic Bethe Ansatz should be exact too (see [0] for evidence 
thereof). We will assume that it is true. 

As is well known [0] the roots of A(u,h) — i.e. the poles of the transfer matrix 
T(u,h), see eq. ( |2.23| ) — are given by the solutions of the Asymptotic Bethe Ansatz 
equations, which only depend on the two-particle phase-shift. In the case of the CSM 
it is hiTsgn(ki — kz). Notice that this phase-shift only depends on the ordering of the 
momenta k\ and This is why these models are interpreted to describe an ideal 
gas of particles with statistics that interpolates between bosons ih = 0) and fermions 
{h — 1). In the dynamical model fl2.1|) the particles have charge and spin. Therefore 
we get two coupled sets of 'Nested' Bethe Ansatz equations — for the general case 
p 7^ 2 there are p equations. They have been presented in [0] : 



hL = hi ^2 7rsgn(fc i - kj) + 



M 



2%Ii -ttJ2 sgn(fcj - A c 



a=l 



N 



it sgn(A a - hp) + 2nJ a = 7r^sgn(A a - h). 



(2.29) 
(2.30) 



i=l 



We have reinstated L, the circumference of the circle to get the dimensions correctly. 
Notice how h drops out of these equations due to the fact that the full Hamiltonian is 
scale invariant (a peculiarity of the ^--type potentials). So rather than sending h — > 
we should let h —>■ oo. There are equations defining the {hi} (one for every particle) 
with charge quantum numbers Furthermore we have M equations defining the 
auxiliary momenta {A Q }. M is the number of particles with a spin {. The {J a} 
are their spin quantum numbers. The J's and J's are distinct integers or half-odd 
integers depending on the parity of and M — N . 

Furthermore we should restrict the Hilbert space to states that only carry spin 
excitations, and no charge excitations (these elastic modes unfortunately don't acquire 
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a gap as h — > O) 4 . As in for instance the lD-Hubbard model we accomplish this by 
leaving the charge quantum numbers in their groundstate configuration and only 
exciting the {J a }. Let us therefore first analyze the absolute groundstate which has 
M = y (for N even), so there are twice as many fc's as A's. The J's and J's are 



consecutive and spaced by one unit. Then eq. ( 2.30|) tells us that between every two 



neighboring A's there must be two fc's. In a spin excited state we have M < y and 
by leaving openings in the J-distribution, we can have more than two fc's sit between 
every pair of A's. Notice that this equation doesn't fix the value of the A's, just 
their positions with respect to the fcj's. From eq. (2.29) we learn that when we order 



the {fcj} such that fcj < kj for i < j: fcj f=a h (i — ^rpj = fc 4 °- There is however 
an 0(x) = 0(H) correction through the A's. Whenever a A sits between two fc's 
they will be drawn together by r. This information is contained in 5 fcj. Now, in the 
same way that the constants of motion are contained in the (D(h)-term in A(u,h), 
their eigenvalues are determined by the O(j^) corrections to the fc's. As A(u,h) has 
eigenvalues Y\^ =1 (u — fcj), T(u, K) = 4- In A(u) = T + hTi + . . . must have eigenvalues 
J2i ^Tfc" = J2i u __ k o 1 _i Sk . ■ So for Ti(u) acting on a state characterized by a set {A Q } 
we find its eigenvalue: 



1 x 



ri(«) = E^TT^E^r'^IAj 



n=0 " i=\ 



oo 1 

= E^TlMAJ- (2.31) 



n=0 



U 



We will now label the A Q 's by m a , their positions relative to the fc's, i.e. if A Q sits 
between fc r and fc r +i, then m a = r. We see that the m a have to be at least two units 
apart since there are at least two fc's between A's according to eq. (|2.30|) . Now writing 
|sgn(A Q — fcj) = 9(A a — fcj) — | {6 is the step function), we have: 

N N ( N + l\ n_1 f M M 1 



4 This is why we don't use eq. (2.12) which gives a much more direct expression for the eigenvalues 
of A (it, h); however, one doesn't know a priori whether an eigenvalue belongs to a pure spin- or charge 
excitation, or a mixture of both. This is something which only a physically motivated method can 
tell us. 
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M 



const + 



a=l 
M 



n) Jt- 

. i=l 



const + e n (m a ) (2.32) 



For small n we can evaluate e n (m a ) exactly: 

ei(m a ) = m a 

t2(m a ) = m a (m a - N) 

, A^ 2 — 1 
e 3 (^a) = 2m a (m a - N)(2m a - N) ^ — m a 

( N 2 -l\ 
e 4 (m Q ) = e 2 {m a ) I e 2 (m a ) H 1 (2.33) 

These results coincide nicely with the numerical values of [0], when we interpret 
the m a as the rapidities of the ISE model! We notice that it is consistent to interpret 
the momentum term P in ( |2.28| ) in Hi and H 3 as J2a m ai i- e - the degree of a polyno- 



mial YHWS wavefunction. Non-YHWS, non polynomial wavefunctions in a Yangian 
multiplet have the same value of P since P commutes with the Yangian Algebra. 

In conclusion we have outlined a method to obtain the constants of motion of 
the ISE model as a strong coupling limit of the CSM with particles with internal 
degrees of freedom. Although the task to actually obtain the invariants is quite 
cumbersome, it can be done in principle. Given the relatively simple structure of the 
invariants we expect there to be some technique that could simplify the computation 
considerably. The construction of integrals of motion presented in this letter, provides 
us with extensive operators that commute with each other and the Yangian symmetry 
algebra. By computing eigenvalues of the invariants through the Nested Bethe Ansatz, 
and comparing them with previous numerical results [0] we provided evidence for 
the validity of the approach. It would be interesting to analyze the cause of the 
miraculous absence of terms containing mixtures of permutations and derivatives. 
Recently a generalization of this method of obtaining the constants of the motion has 
been generalized to spin models with a trigonometric rather than a rational i?-matrix 
in the FCR eq. ( 1.18Q for their transfer matrix [0]. These trigonometric analogs of 



Hise do not have the Yangian Yfafa) but rather the quantum group U q {sl2) as their 
symmetry algebra. 
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Chapter 3 



Dynamical T = correlations of 
the ISE-model in a magnetic field 

3.1 Introduction 

In this chapter we will devote ourselves to computing an observable in the ISE-model: 
the dynamical spin-spin correlation function or Dynamical Structure Functions (DSF) 
in a non-zero external magnetic field. If the ISE-model were to be realized directly 
in the nature, the DSF would be of foremost interest since it measures effectively 
the intensity for inelastic neutron scattering off the spin-chain. A more important 
reason to be concerned with it is that, because of the simplicity of the model, we 
can actually obtain exact expressions for the DSF. As a matter of fact, the zero 
field case for this model was the first strongly interacting system for which such 
correlators could be obtained at all [0, 0]. The static structure factor had already 
been obtained in 1988 [0, 0], based on work by Sutherland from the 70's [0]. It turns 
out that both have all the qualitative features of the complicated NNE mode, but are 
quantitatively much cleaner. The underlying reason for this cleaner structure is that, 
where normally a neutron, colliding inelastically with the ring of spins, tends to excite 
many quasiparticles (although high numbers of them are less likely), in the ISE-model 
there can be no more than four, and in zero field just two. This among other things, 
enables the aforementioned exact approach. In this chapter we want to generalize this 
to non-zero magnetic fields. The rationale is that at T = and B = the groundstate 
contains no elementary excitations, the spinons; by tuning the magnetic field, which 
acts as a chemical potential for the up-spinons, we can select any non-zero number of 
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them in the groundstate. This allows for a richer set of excitations and also enables 
us to make contact with the excitations over the ferromagnetic state (in high fields) 
which should be magnons — with spin | rather than 1 [0]. 

In the next section we will define the dynamical structure function and its proper- 
ties in the ISE-model. In section we will discuss the various contributions to the 
DSFs. Finally, section |3.4| compares our results to those of Miiller et al. [0]. They 
obtained approximate analytical expressions and numerical results for dynamical cor- 
relators in the NNE chain. 



3.2 Dynamical Structure Functions 

The DSF for a spin model at zero temperature is defined as: 

(GS|5^(*)S£(0|GS) = S ij (m,n,t-t'). (3.1) 

|GS) represents the groundstate of the system, and S m (t) is the local spin operator, 
acting at site m and time t. Its time-dependence is implemented via the Heisenberg 
representation: S m (t) = e iHt S m e- iHt . The DSF measures the response of the system 
to flipping/imposing a certain spin on a site n in the groundstate at time t' and 
measuring the effect on a site m at time t. Since H conserves S z we cannot create 
an up-spin somewhere and have it evolve into a down-spin. In other words, there are 
restrictions on the possible choices of i and j in (|3.lQ . We find that only S~ + , S + ~ 
and S zz are non-zero {S xx = S™ = \ [S+~ + S~ + }, S 3 * = ± [S~ + - S+-}). 

The DSF fl3.1|) tells us what kinds of excitations a local spin operator can generate 
in this system. This is obvious, once we insert a complete set of eigenstates {|z^)} 
between the two operators: 

5 ab (m,n;t-t') =E( GS l 5 mk)(^l^|GS)e-s ( *'- i)( ^- Eo) . (3.2) 

V 

This can be simplified further in the case that the model is translationally invariant 
(such is the case for both NNE and ISE models), i.e. both |GS) and \ v) are momentum 
eigenstates. In that case S ab will only depend oam-n and we can consider its time 
and space Fourier transform: 

S{Q,E) = J2KS(E-(E„-E Q ))5(Q-(p u - Po )) 

M a v = 27r|(z/|S a (Q)|GS>| 2 , a = ±,z, (3.3) 
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and S(Q) = ^= J2n=i e_m<9 S n , and (p v ,E v ) and (p , E ) are the momentum and 
energy of the state |z/) and |GS) respectively. Notice that the support S(Q, E) is zero 
except when (Q, E) corresponds to the excitation energy and momentum of a state 
contained in S(Q)|GS). 

For a general model, this set is constrained by global symmetries, such as mo- 
mentum conservation, and the Wigner-Eckart theorem — from the vector character of 
S(Q). For the ISE-model, however, the number of excited states contributing to the 
sum (|3.3| ) is very limited, which is what we alluded to in the the introduction, when 
we mentioned the simplicity of this model. We will find that the relevant excited 
states differ from the groundstate only in the following way: in low fields the number 
of spinons can be the same as that in the groundstate, or it can increase by two; 
in high magnetic fields, there are two additional excitations: a right or left moving 
magnon, or the number of spinons can decrease by two. 

The reason for the different low and high field behavior is that the character of the 
groundstate changes with magnetic field. This is caused by the additional Zeeman 
term in the Hamiltonian: 

Hzeeman = ~h^2 S- (3.4) 

i 

At h = the groundstate lies in the M = y sector (M is the number of down-spins). 

As we increase the field, there comes a point where h is equal to the h = energy 

difference between the (absolute) groundstate in the M = y sector and the one in 

the M = y — 1 sector, and the energies of these states cross. We find that for any 

value of h the groundstate is always YHWS. 
This is easy to see, if we consider the expression 

for the energy of a Yangian multiplet in terms 
of the pseudomomenta ( |1.45| ), E oc J2i e ( m i) — 
\Hi m i( m i ~ N). We see that adding pseudo- 
momenta lowers the energy. 

I.e., the Yangian multiplet with one less 
pseudomomentum, and S z higher by +1, has 
higher energy. So in a given sector of S z the 
groundstate is a member of the subset of YHWS (with that value of S z ). From 
the figure it is clear that in the groundstate the pseudomomenta will arrange them- 
selves symmetrically around m = y. The expression for this state is actually known 












N 








3» 

m 
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analytically [0]: 

M 2 M N 

il>G S (n 1 ,...,n M ) = II { z ni ~ z r h ) 11^ + ( 3 - 5 ) 

i<j i=l 

Although the groundstate magnetization S z = \{N — 2M) is a discontinuous function 
of h, it becomes smooth in the thermodynamic limit. With the set of pseudomomenta 
{ mi = K-M + 2i-l}fi 1 , the groundstate energy is E (M) = M(4M2 ~ 4 3jv2 ~ 4) - 

§(-/V — 2M). Minimizing with respect to M gives: 

As expected, with the magnetic field, the number of fully polarized spinons N sp = 2S Z 
in the groundstate grows. For h > h c the groundstate is completely ferromagnetic. 

For the sequel it is useful to have expression of energy and momentum of an 
eigenstate of Hi SE in terms of "spinon momenta" as well [0]. If we label the M + 1 
orbitals from right to left by spinon momenta — k < k < k = spaced by ^ we 
get: 

P = k + Nko mod 2n, 

k 

E = Y / e(k)n k( r + ^ E V(k-k')+ E(M,N), (3.7) 

ku kk',aa' 

where e(k) = —(k^ — k 2 ), V(k) = v s (ko — \k\) and n^ c is the number of spinons in the 
orbital with momentum k and spin a [0]. E(M, N) only depends on the total number 
of sites and spinons. We also recognize v s as the spinon velocity at the zone 
boundary, ko, in the groundstate. Notice from the spinon and magnon dispersion 
relation, we see that the "mass" of the magnon is (-2) times that of the spinon. The 
system is gapless, so this mass doesn't refer to the energy that is required to create 
a pair of them, but to the curvature in their band. 

Before we turn our attention to the precise selection rules that restrict the set 
of accessible intermediate states, let us first discuss an intuitively more attractive 
(overcomplete) basis of the Hilbert space: the set of localized spinon wavefunctions 
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[0]: 

M M N sp 



y(ni,...,n M \a h ...,a Nap ) = J[ (z ni - z n .) \{ z ni \{ \{ [z ni - z a .) . (3.8) 

i<j i=l i=l j=l 

The {on} are the locations of localized spinons that can point up or down; the {n;} are 
the positions of down-spins (other than those of possible localized spinons pointing 
down). Notice that the wavefunction prevents the {n;} from coinciding with the 
spinon sites. We call the complement of the set of spinon sites the condensate. It is a 
singlet under the action of total spin. Furthermore for N sp = eq. (|3.8|) represents the 
exact groundstate wavefunction for h = 0. The usefulness of these states is limited 
by the fact that they are not mutually orthogonal and worse, overcomplete. However, 
based on numerical evidence for up to 12 spinons, it is clear that the space spanned 
by states Q3.8| ) with a fixed number N sp of localized spinons contains only eigenstates 
of the Hamiltonian belonging to Yangian multiplets with N sp or less spinons. 

The subspace of states that have a fixed number of N sp localized spinons, all po- 
larized, has the pleasant property that it only contains eigenstates of H with precisely 
N sp spinons. This is clear from the fact that these localized spinon wavefunctions — 
although not eigenstates of the Hamiltonian — are annihilated by both Jq and Ji (see 
Appendix A). This is consistent with the fact that fully polarized spinon eigenstates 
are supposed to be of a polynomial form with degree < A in the {z ni } [0], just like 
wavefunctions (|3.8| ). In that same reference it is also shown that these fully polarized 
localized spinon states are complete as well and span all YHWS. 



3.3 Selection Rules for Spin operators. 

We continue discussing the various Yangian multiplets in terms of their pseudo- 
momenta — the zeroes of /2(w) — rather than the Drinfel'd zeroes. It is customary 
to denote the set of pseudomomenta by a sequence of A — 1 O's and l's, where a 1 



indicates a pseudomomentum (the locations of the big crosses in fig. |1.1| ), and O's (the 
small crosses and circles) pad the l's so that all integers in the range are covered. 
Such a binary sequence is also called an occupation pattern. 

We will now discuss the relevant selection rules for the matrix elements M% in eq. 
( |3.3|) . From numerical evidence up to N = 16 sites it has become clear that there are 



Chapter 3: Dynamical T = correlations of the ISE-model in a magnetic Geld 62 



only an unexpectedly small number of those matrix elements that are non-zero. To 
resolve parity and other accidental degeneracies between the Yangian multiplets we 
split these degeneracies 1 by actually diagonalizing H + XH 3 , where H 3 is the second 
integral of the motion for this model as presented in eq. ( |1.14| ). The eigenvalues 



of this operator allowed the Yangian occupations sequences to be unambiguously 
identified. States in the multiplets are partially resolved by fixing S z and 5**°' (a 
unique resolution of states would be obtained by adding another term /iJ • Ji oc Jq 
to the Hamiltonian 2 . This would correspond to a basis of states within the Yangian 
multiplet obtained through the Algebraic Bethe Ansatz [0]). 

Now let us denote the eigenstates of this model as |T, /z) where T labels a Yangian 
multiplet through an occupation sequence and fi labels the state within the multiplet. 
Furthermore define M FF ' (S^) ^ = (T/j,\ \Tfjf). Then the observation made above 
implies that the matrix M rr (Sf) vanishes if the occupation sequences T and V differ 
"too much" in a sense made precise below. This situation is analogous to an ideal 
gas, where if O is a one body operator: (a\0\/3) = if the occupation number 
configurations of |a) and \(3) differ on more than one orbital. E.g. O = p(x) = 
J2kk' e l ^ k ~ kn>a c\.Ck> can add or take out a single particle from an orbital in an ideal gas, 
but in an interacting gas it could add unlimited numbers of particle/hole pairs. 

We have formulated a set of empirical rules, summarizing numerical studies on 
small systems with iV < 16. We expect that these will be derived when a complete 
understanding of the model is developed. According to these rules, if two multiplets 
T and r" do not differ too much, there will always be a pair of states fi and n! in 
either multiplet for which the spin matrix element is nonzero. 

EMPERICAL RULE 1: If ir(T,m,n) is the total number of ones in positions m 
through n in Yangian occupation sequence T then: M rr (S?) 7^ Q iff. |7r(r,m, n) — 
7r(r', m, n)\ < 1 for any \ <m<n<N— 1. 



The rule is illustrated in Fig. |3.1|. A general consequence of this rule is that when we 



choose m = 1 and n — N — 1, it follows that the total number of ones in a sequence 
can't change by more than one, i.e. the total number of spinons can only change 



1 Due to the fact that the expression of the energy of a state (1.45) is sum of squares of integers. 
2 This term splits the spin multiplets within a Yangian occupation pattern, with identical Srbt 



but different rapidities {A^} in eq. (1.26) 
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(r = ioooiooioooo, n\ s? \r = 100010001000, //) ^ o 
(r = ioo oiooi oooo, n\ s? |r" = 001000001000,/i") = 

Figure 3.1: According to the selection rule, there will be states \x in multiplet T that 
are connected to others jJ in V through a local spin operator. In V and T" there are 
none, e.g. since |tt(T, 4, 8) - 7r(F, 4, 8) | = 2 - = 2 > 1. 

by +2,0, or -2. It is remarkable that according to the rule this also holds on any 
corresponding subsequences of the occupation number sequences. 

The zero magnetic field DSF has been computed in [0]. The particular structure 
function computed happened to be S~ + (Q,E) (the others are identical because of 
rotational invariance). This function is governed by excitations present in |GS). 
Since the zero field groundstate contains no spinons the rule tells us that we can 
only expect excitations with or 2 spinons. As |GS) has S z = +1 the former is 
ruled out and in the states in the multiplets with 2 spinons, both must be polarized. 
This was to be expected since we can expand |GS) in a set of localized spinon 
wavefunctions ( |3.8| ) containing two polarized spinons: 

n— l 2 1 — (— ) m 

(n 1 ,...,n M - l \Sf |GS) = £ — - — K -^-V(n 1 ,...,n M - 1 \i,i + m). (3.9) 

m=l N 1 — Z m 

So one of the spinons seems to be sitting on the site on which acted and the other 
is an even number of sites removed from it. 

Since we know the spinon dispersion relation, we can demarcate the support of 
S~ + (Q,E) in the (Q, _E)-plane in Fig. |3~2"| . The main steps of the computation in- 
volving the matrix elements M + , i.e. the weight of the DSF at a point (Q, E) on the 
plot, are the following: since |GS) only contains states with two fully polarized 
spinons, it must be built out of YHWS. These wavefunctions are functionally identi- 
cal to eigenf unctions of the Calogero-Sutherland model at coupling A = 2 of particles 
moving on a ring. Since both wavefunctions are of a polynomial form with degree < iV 
the computation involving a sum over sites is identical to taking an integral over the 
ring in the continuum model: Jjf dne lkn = J2n=o e%kn $ k & { — (N — 1), . . . , N — 1}. 
The action of in the spin chain is translated into a particle destruction operator 
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Figure 3.2: The shaded region show where S h (Q, E) is nonzero for h = 0. The top 
boundary corresponds to excitations with two spinons that have identical momentum; 
on the bottom boundary one of the spinons has fixed momentum ±7T. E is given in 
units of (f ) 2 ^. 

ty(x,t) (Si~ removes a down spin). So the S~ + (Q,E) DSF reduces to the Greens 
function (GS| ^(x, t)\I/(0, 0) |GS) in the Calogero-Sutherland model. It can be com- 
puted in the thermodynamic limit, in which case it can be mapped unto a Gaussian 
hermitian matrix model correlator. The result is: 



by vi,V2, which are quickly identified with the velocities of the two spinons in the 
excited state. The (ty^) Greens function has been obtained recently at finite N as 
well: both the groundstate, the YHWS \u) it connects with, and the action of Sf are 
expressible as symmetric Jack-polynomials, and the matrix element is essentially the 
Clebsch-Gordan coefficient of the three [0]. 

For /i / the three different structure functions 5* h (Q,E), S ZZ (Q,E) and 
S + ~(Q,E) will not be equal, since |GS) is no longer a singlet. In fact it has 
S z = s tot = = where Ngp = N^h) is given by (Q. This difference be- 
tween the three correlation function is also expressed in two additional global SU (2) 
selection rules [0], which rule out certain matrix elements based on the total spin S tot 




(3.10) 



with Q = —-^-{vi+V2) and E = e(vi)+e(v2)- The DSF matrix element is parametrized 
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and S z of the final state \i inside the Yangian multiplet V. 

In the first place there is the Wigner-Eckart theorem for vector operators such 
as Sf, which tells us that in order for (T, Sf |GS) to be nonzero the total spin S 
of \T,fj,) must satisfy S — 1 < S < S + 1. Secondly, for any state in a multiplet 
with N sp spinons, we have > S > S z — where equality only holds for the YHWS, 
which has all its spinons up. Now put this together with the fact that Sf raises S z 
by +1, lowers it by 1, and S z leaves it the same. Classifying states according to 
their S z and S tot as types (i) - (vi), following Miiller et al. [0] we find the following 
contributions: 

Sf | GS ) contains states with S = S + 1 and AA^ sp = +2 (type (in) ). 

S z |GS) contains states with S = S + 1 and AN sp = +2 (type (i) ) or S = S and 
AN sp = +2,0 (type (ii) a > b ). 

S^ |GS) contains states with S = S + 1 and AN sp = +2 (type (iv) ) or S — Sq and 
AN sp = +2, (type (v) a ' b ),orS = S Q -l and AN sp = +2, 0, -2 (type (vi) a ' b < c ). 



Since we have an additional quantum number to label states: N sp , we added Latin 
superscripts a, b, c to the Roman numerals. All 10 contributions are summarized in 
table (HOD. 

We will now investigate all three structure functions individually following these 
selection rules. 

3.3.1 S (Q,E) 

Type (in): AN sp = +2, AS = +1. 

For the occupation sequence of the groundstate in a given magnetic field (e.g. 
000010101010000) let us label the zeroes in the leftmost orbital as the left spinon con- 
densate and the ones in the rightmost orbital as the right spinon condensate. From 



table |3.1| we learn that action of Sf on the groundstate only produces states with 
two more spinons, i.e. one less 1. This 1 has to come out of the center • ■ • 10101 • ■ ■ 
region. We can't take more than a single 1 out of the center region — and place 
it into the left or right spinon condensate — since this would imply a violation of 
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Table 3.1: List of matrix elements contributing to the DSF. The correlation functions 
in the last column refer to those in the Calogero- Sutherland model as they appear 
in eqs. (|3~T0| , gl3j gll [3~17| , |3~20| , |3~2TD . The entries marked with a (f) do not 



survive in the thermodynamic limit. For the example occupation sequences we act on 
the specific groundstate: 000|l0101010l|000, with the symbols | delimiting the center 
region. In the column under "Excitation" S denotes a spinon and M. denotes a 
magnon. 



DSF 


Type 


S tot 




AN sp 


Typical contributing 
Yangian multiplet 


Excitation 


Matrix element 


S~+(Q,E) 


(Hi) 


s = s + 


1 


AN sp = 


+2 


000|100100101|000 


25 




(5'=5o + l) 






















S"(Q,E) 


(*)t 


s = s + 


1 


AN ap = 


+2 


000|100101001|000 


25 




(S'=S D ) 


(iiy 


s = s 




AN ap = 


+2 


ooojiooioiooijooo 


25 




(uf 


S = S 




AN sp = 


- 


100|101001001|000 


25 + M 


(pp) 


S+-(Q,E) 


(iv)^ 


s = s + 


1 


AN sp = 


+2 


000|100100101|000 


25 


(5o+1) J 250+1) <* t *) 


(S'=S -1) 


«t 


S = S 




AN ap = 


+2 


000 


010010101 


000 


25 






(vi) a 


S = S - 


1 


AN ap = 


+2 


000 


010010101 


000 


25 






(vf 


S = S 




AN ap = 


= 


000 


100100101 


010 


25 + M 






(vi) b 


s = s - 


1 


AN sp = 


= 


100 


100101001 


000 


2S + M 




(vif 


S = So - 


1 


AN sp = 


-2 


100 


100100101 


010 


2S + 2M, M 
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Rule applied to the center region. Taking out a 1 in the center region is equiv- 
alent to inserting 2 spinons there. A typical nonzero matrix element would be 
(0001001001000, fi\ St |0001010101000), with the two spinons residing in orbital two 
and three. 

All this means is that we get a simple two spinon spectrum, as in the zero magnetic 
field case. The only difference is that now the momenta of the spinons can only vary 
from — ko to k , where k = jjj(N — N sp ) decreases with increasing magnetic field 
as N sp = N sp (h) according to eq. (|3.6j ). The support of S~ + (Q,E) is essentially a 
squeezed version of Fig. EJ72, 



As for the weight associated with 2-spinon excitations: the calculation for the 
zero magnetic field case (|3.10| ) carries over without problems. The reason for this is 
twofold. In the first place, the nonzero magnetic field groundstate wavefunction (|3T 



is of the same Jastrow-form as the zero field one, with just an extra phase factor 
fli appended. When we take the matrix element, the phase factors from ket 

and bra part cancel each other. Secondly, the excited states are again of the YHWS 
type (S = ^p) and have to be polynomials. The mapping onto a Calogero-Sutherland 
model matrix element remains therefore legitimate. 

The contribution of just two-spinon YHWS excitations was to be expected since 
we can expand any fully polarized localized spinon wavefunction with N sp spinons 
acted upon with in terms of a set containing N sp + 2 spinons: 



(^ ai )(nir--,%-i) = 

II _ \ ^( n ii ■ ■ ■ ,n M -i\ai, . . .,a Nsp ,i,p). (3.11) 

Here V is a set of M random sites on the circle excluding the spinon sites {a{\ (see 
Appendix B). Eq. ( |3.9| ) is a special case of this expansion with V equal to the sites 
that are an even number of steps removed from the site on which the local spin 
operator acts. We could also have realized that the number of spinons can't go up by 
more than 2 when we consider that S + and annihilate |GS), indicating that 
the latter must consist of purely YHWS with S z = So + 1! 
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3.3.2 S ZZ (Q,E) 

For the DSF S ZZ (Q,E) we find similar simple excitations that contribute, although 
at present not all resulting matrix elements can be computed. From the combined 



selection rules we find three types of excitations: 



in) 



and 



n 



in table |3.1| . 



They all have in common that AN sp = or +2. This isn't surprising: let us consider 
the state Jf (Sf |GS )) . This state is annihilated by Jf and S + , so it must be YHWS 
(see Appendix C). Therefore, since the first action of Ji doesn't change the number 
of spinons, S? |GS) must be a mixture of states that contain not more than 2Sq + 2 
spinons, where 2 So is the number of spinons in the groundstate. 

We will now discuss the individual types and where possible compute the values 
of the matrix elements. 



Type 



AN. 



sp 



+2, AS 



-1. 



Having identical selection rules, these states sit in the exact same Yangian- and 
spin multiplets as the type (Hi) states, and therefore they also contain two additional 
spinons. However their S z = S tot — 1, so they are no longer YHWS like the type (Hi) 
states. Nevertheless they are related by a simple application of S~: 



IrA/v -+2 s * ot=s o+i 

I J. ,ZY/V SJJ — s z =S 



) 



'—S" ir), 

Sn+l) 1 " 



V 2 ( 5 + 1 ) 



(3.12) 



where |T) denotes the YHWS of the multiplet with occupation sequence T of type 
(Hi). This allows us to reduce a type (i) matrix element to one that is a multiple of 
a type (in) given in eq. (|3.10| ): 



(r,Aiv sp = 



=So+i 



| Sj | GS ) 



i 

2(S + 1) 

i 

2(50 + 1) 



r|[s + ,s 4 2 ]|Gs; 



(r|s+|GS) 



(3.13) 



where we used that S + annihilates the groundstate. 



Type 



ii) 



AN, 



sp 



-2, AS = 0. 



This type of state is a member of the same kind of Yangian multiplet as type (Hi) 
and (i) — i.e. with two extra spinons — but it sits in a spin multiplet that does not 
contain the YHWS. States of this kind are supplied by the Algebraic Bethe Ansatz, 
eq. (|L2§): T-(A )|r). A is given by eq. ( |Q7|) for n = 1: P{ p X ^ ] = 1, where P(u) 
is the Drinfel'd polynomial of occupation sequence V. This equation has a number of 
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Figure 3.3: Graphical illustration of the spin multiplets accessed in ^-excitations of 
types (i) and (ii) a . 

trivial solutions, when Ao corresponds to two adjacent Drinfel'd roots. These trivial 
solutions are to be discarded (they correspond to an orbital with more than one spinon 
in it; these spinons can only be in a symmetric state, and there is a reduction in the 
naively expected number 2(5*0 + 1) — 1 of £ = So states). The problem is, how to 
commute the T~(Aq) through Sf in the matrix element (GS|£fT~(Ao)|r). 



Type (ii) b : AN, 



sp 



0, AS 



0. 



Since this type of state has AiV S p = the number of ones in its occupation sequence 
must be identical to that in the groundstate. As with the previous three types we can 
delete just a single 1 from the center, leaving behind two spinons pointing up. This 
1 then must go into either the left or the right spinon condensate. Therefore a typi- 
cal nonzero matrix element would be: (00100|l00100l|00000| Sf |00000|l01010l|00000) 
where the | just helps to draw attention to the center region. Rule [I] rules out any 
additional ones leaving the center region. The additional 1 on the left or right can be 
interpreted as a magnon with S 2 = — 1. The limiting case where the magnon "fuses" 
with two spinons at a boundary between a condensate and the center region gives us 
the groundstate. 

These states must be YHWS since they have = S = S z , like |GS). This fact 
allows us to calculate the corresponding matrix elements. Since both groundstate and 
excited state are YHWS, a mapping onto the Calogero-Sutherland model is valid. In 
this case we need a groundstate density- density correlator (GS| p(x, t)p(x f , t') |GS) 
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since S z measures the presence or absence of a down-spin (i.e. a particle in the CS- 
model). This calculation has been done by Altshuler et al. [0], in the thermodynamic 
limit, by studying the repulsion of energy levels in a random matrix model under 
a varying perturbation. The energy levels are identified with the positions of the 
particles and the strength of the perturbation corresponds to imaginary time. Their 
original expression depends on three parameters (called A, Ai and A 2 ), the latter two of 
which are compact, and the first one is unbounded. This is precisely what we expect 
from our selection rule: 2 spinons restricted to the center region with momenta in the 
range — k . . . + k and a magnon that can go off all the way to the right or left (i.e. 
±oo in the thermodynamic limit). In terms of the velocity, v, of the magnon — with 
dispersion relation ( |1.45| ) — and spinon velocities vi,v 2 their result is as follows: 



\(v,v 1 ,v 2 \ p{Q) |GS)| 2 oc 

(u-t; 8 )(i; + t; 8 )|i;i -v 2 \ 

^J(v 2 s -vl)(v 2 s -vl){v - v^iv - v 2 y 



[{v - v x ) + (v - v 2 )f . (3.14) 



E = e(t>i) + e(v 2 ) + e(v) and Q = ^-(v — |(t>i + v 2 )). For completeness we give the 
relation between the v 's and the A's : 



A = v/v s 



(1-A?)(1-A 2 ) = (^) 2 . (3.15) 

Numerical data for S ZZ (Q, E) can be found in fig. (|3.4| ) at values of h close to and 
h c . Fig. ( |3.5| ) shows the corresponding support of S ZZ (Q,E) in the (Q,E) plane, as 
predicted from the selection rule, in the thermodynamic limit. We notice that for 
finite-size systems some of the features near the lower boundary are absent. 



3.3.3 S+ (Q,E) 

Finally we discuss the S + ~(Q, E) DSF which is governed by the excitations of types 
(iv) - {vi) c that are present in S~[ |GS). 



Type (iv): AN sp = +2, AS = +1. 
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This type of state is very similar to types (Hi) and (i), as a matter of fact they 
all reside in identical Yangian- and spin multiplets. Therefore type (iv) states differ 
from the groundstate only by two extra spinons in the center region. They are related 
to their YHWS |T) (of type (Hi)) through: 

|r,A^ +2 , -Eg; > = 7 — k — (S- f |D. (3.16) 

Analogous to the calculation for type (i) states, a matrix element of type (iv) can 
now easily be reduced to one involving type (Hi): 



(r,AAf sp =+2, 



S™ Z =S +1 I qz ipql 

s z =s -i \°i I 00 / 



(45 () +2)(25 () +2) 



(2S + 1)(S + 1) 



(T\(S + yS- |GS> 



T|5+|GS) 



(3.17) 



The last matrix element in this equation has already been computed for the S h (Q, E) 
DSF. However the energy of the corresponding excited state here is shifted by 2h in 
comparison, because of the Zeeman term in the Hamiltonian. 



Types (v) a and (vi) a :AN sp = +2, AS = 0, -1. 
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Figure 3.5: The regions where S ZZ (Q,E) is nonzero for low and high h. The area 
shaded dark contains the contributions from the excitations with AN sp = +2 (types 
((i) and (ii) a ). The excitations with AN sp = — type (ii) b — can also occupy the 
lightly shaded area. 



States of types (v) a &nd(vi) a contain a two-spinon excitation as well — like types 
(Hi), (i) and (iv). However, since they don't reside in the spin multiplet of their 



YHWS 



-N. 



> S), the associated matrix elements are unknown. Type (v) can be 



related to those of type (ii) a , though. 

2 



, qZOZ q 

(r,AiV sp =o, 5 s z=So s ° 



\ I^Igs) 



i 

2S 



;r,Aiv sp =o, s s r = =/ ° |5?|GS> 



(3.18) 



Type (v i) b is even worse, since it involves the application of two T (A) operators. 
Types (v) b : AN sp = 0, AS = 0. 



Type (v) b states are very similar to those of type (ii) b , they only differ in S z by 
-1. Therefore both contain two excited spinons and a single left- or right moving 
magnon. They are related to each other by: 



|r,Aiv sp =o, s , = -_° ) 



(3.19) 



- v — u 

and |T) is the type (ii) b YHWS of the multiplet with occupation sequence T. We can 
— " trivially relate the matrix elements of type (ii) b and (v) b : 



now 



s z =s n 



1 

2S„ 



(T| [S+,Sr] |GS) 



^|<r|5?|GS>| : 



(3.20) 
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The last matrix element is listed in eq. ( 3.14Q . 



Type (vi) b : AN sp = 0, AS = -1. 



These states reside in the same Yangian multiplets as the previous type however 
they are not in the spin multiplets of the YHWS: -| £ > S = S — 1. The problem is 
very similar to that for case (ii) a . 

Type (vi) c : AN sp = -2, AS = -1. 

This last class of states has AN sp = —2 and therefore the number of ones in its 
occupation sequence goes up by one compared to the groundstate. The selection Rule 
[l] only allows the extra 1 to go into the left or right spinon condensate. As before we 
are also allowed to take a 1 out of the center region and bring it into the left or right 
spinon condensate. Notice however that the rule forbids both of the ones to go into the 
same condensate: one has to be left moving and the other must be right moving. The 
result is an excited state with two magnons and two spinons. A typical nonzero matrix 
element would be (0100|l0010100l|0001 S[~ |0000|l0101010l|0000). The | just helps to 
guide the eye. Also present are YHWS of multiplets from the limiting cases where one 
of the magnons fuses with the two spinons in the center; this leaves a multiplet with 
nothing but one magnon; example: (0000|l0101010l|001o| S^ |0000|l0101010l|0000). 
This single magnon excitation is familiar from the strong field regime. 

Since type (vi) c states are YHWS C^f = S), as is the groundstate, we can 
repeat the calculation of the matrix elements by a mapping onto the Calogero- 
Sutherland model. Because S~ creates a down spin, it corresponds to a particle 
creation operator in the CS-model. The relevant correlation function is therefore 
(GS| *&(x, t)^(x',t') |GS). As in the (tfttf) case for type (m) states, a further map- 
ping onto a Gaussian Hermitean matrix model allows one to calculate the Fourier 
transform of this correlation function [0]. The result consists of two parts and 
q(2)_ q(2) j g p aram etrized by four variables, two of which are compact: \vi\, < v s , 
and two are noncompact:t> < — v s , v' > v s : 



[v,v',v 1 ,v 2 \S (Q)\GS) 



2_ J (V - Vi)(v' - V 2 ) + (V - Vi)(v - 

\ (v - V ') 

Q 2 - V 2 s ){v' 2 - V 2 s )\Vl - V2\ 




(3.21) 
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Energy and momentum in terms of the u's are given by Q = -jj-{v+v' — 2(^1 + ^2)) an d 

Vg). It is obvious that the compact parameters are to 



be identified with the spinon velocities and the non-compact ones with the magnon 
velocities. is parametrized by just one variable \v\ > v s , which is noncompact 



HS-(Q) |GS> 



1 V - Vs 

2n v + v s 



(3.22) 



This contribution corresponds to the case where one of the magnons has \v\ = v s , 
and "fuses" with the two spinons, thus leaving but one magnon. Since we now know 
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Figure 3.6: Non-zero region of 4S XX {Q,E) = S + -(Q,E) + S- + (Q,E) in the 
(Q, £7)-plane for low and high /i. The contributions of AiV sp = +2 (types (if), (t>) a 
and (vi) a ) live in the dark shaded region. These excitations will survive in the limit 
h — > 0. The AN sp = 0, —2 can also occupy the area shaded light. For h — > h c only the 
1 magnon contributions survive, as can be seen in Fig. |3.7| . These high field magnon 
excitations are indicated by the thick lines. 



all possible excitations contributing to S + (Q, E) we can draw its support in Fig. ^]6| 
for low and high values of h. Fig. |3.7| 3 show numerical data on S + ~(Q, E) for those 



values of h. Table ^3] summarizes the selection rules and available information on 
matrix elements. 



I would like to thank Dr. K. Lefmann and Dr. C. Rischel for pointing out an error in this figure, 
in the original publication. 
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S (Q,E)atlowh 



S (Q.E)athighh 




Figure 3.7: S XX (Q, E) for low and high h, from numerical diagonalizations on a 10-site 



chain. E is in units of (4£ 



\2v s 



3.4 Comparison to the Bethe Ansatz Model 



In 1980 Miiller et al. [0] did a similar calculation of DSFs for the nearest neighbor 
Heisenberg chain. They identified certain types of states called spin wave continuum 
states (SWC) as carrying the dominant contribution to the DSFs. These SWC states 
can be described in a Bethe Ansatz rapidity language by occupation sequences, just 
like the Yangian multiplets in the ISE model. As it turns out, these SWC states 
correspond to exactly the same rapidity sequences that are favored by our selection 
rules in the ISE-model! Although the dispersion relations for the NNE Bethe Ansatz 
rapidities are different from those in the ISE-model, the authors find the support of 
the DSFs in the nearest neighbor model to have essentially the same shape as we do 



in the the inverse exchange case, see fig. |3.8| . 

However, in the nearest-neighbor Heisenberg chain there are some "anomalous" 
states, characterized by a change of more than ±2 spinons, which contribute to a 
lesser degree to the structure functions and don't lie within the bounds found by 
the authors. In the ISE-model these contributions are completely absent and once 
again we find this model to have a surprisingly clean structure. So in this sense the 
ISE-model is an ideal spinon gas, whereas in the NNE Heisenberg chain the spinons 
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interact. 



See figure 15 in Phys. Rev. B 24, 1429 (1981) 



Figure 3.8: The S zz structure functions according to Miiller et al., Phys. Rev. B 
24, 1429 (1981), for iV = 10 in low (S^ s = 2) and high (S^ s = 4) magnetic fields. 
Notice the good match with its ISE counterpart fig. [3.4| , except for the low weight 
NNE states outside the 2-spinon/l-magnon region. The solid lines are boundaries of 
regions that receive contributions from just 2-spinon excitations. Their dispersion is 
trigonometric, rather than quadratic, but qualitatively similar. 

The characteristic shape of the region in the Q, E plane that carries all (ISE) or 
most (NNE-model) weight has recently been seen experimentally too [0]. Tennant et 
al. scanned along a line in the energy-momentum region with the help of neutrons 
scattering off KCUF3 at low (but not too low) temperatures, see fig. |3.9| . The intensity 
obviously peaks, when traversing the "allowed region". In the forbidden regions the 
intensity drops, but not all the way to zero. This might happen because KCuF 3 is 
better described by a NNE model than by the ISE-model, but the fact that T 7^ can 
cause a more smeared distribution too. Finally, there may be non-magnetic scattering 
centers that contribute to the structure function as well. 
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See figure 1 in 

Phys. Rev. Lett. 70, 4003 (1993) 



Figure 3.9: Neutron scattering data from a quasi one dimensional spin chain KCuF 3 
according to Tennant et al. Phys. Rev. Lett. 70, 4003 (1993). Response is large in 
regions where the incoming neutrons can excite just two spinons. 



Miiller et al. also gave general rules, based on comparing Clebsch-Gordon coeffi- 
cients, determining which matrix elements will survive in the thermodynamic limit. 
Their conclusion is that the only surviving ones have Srot = S z . This means only 
excitations (Hi) for S' + (Q,E) , (ii) a and (iif for S ZZ (Q,E) and (vi) a , (vi) b and (vi) c 
for S + ~(Q,E) remain relevant. (Exceptions are single excitations with Q = 0, since 
these correspond to S + , S z and S~ which give macroscopic contributions, as we can 
see in the figures). These contributions are trivial to compute. 



3.5 Conclusion 

We found a remarkably simple selection rule for nonzero matrix elements of local spin 
operators between eigenstates of the ISE-model, which is reminiscent of the ideal gas 
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single-particle selection rules. One of the consequences of this rule is that the total 
number of spinons can only change by 0, ±2. Within the occupation sequences this 
holds locally as well. 

In the particular case that one of the states in the matrix element is also the 
groundstate in a magnetic field (i.e. fully polarized spinons, condensed into the left- 
and rightmost orbitals in equal amounts), the general selection rule only allows ex- 
citations with no more than 2 spinons (the rule applied to the center region) and 
one left- and one right moving magnon (the rule applied to the condensates on the 
left and right) 4 . This implies that the structure functions based on the matrix ele- 
ments involving these states have a finite support in a region dictated by convolving 
the dispersion relations of these particles (Figs.[T5] and |3.6|) . These regions have the 
same shape as those for the nearest neighbor Heisenberg chain. The latter model 
carries some weight outside these regions as well. Therefore the ISE-model has a 
much cleaner spinon structure than the NNE-model. 

Matrix elements that connect a number of the states in these regions to the ground- 
state through the local action of a spin operator have been presented. However, in- 
formation is lacking on those types that involve states not in the spin multiplet of 
the YHWS. This is particularly bothersome for types (ii) a , (vi) a and (vi) b since these 
will survive in the thermodynamic limit. Their calculation would allow a full recon- 
struction of the h 7^ DFS for the ISE-model. A more algebraic treatment involving 
Yangian operators should provide more insight. 

We believe that eventually all dynamical correlation functions of the model will 
be explicitly found, namely (S(Q,E)) at finite magnetic field and temperature. This 
is a goal that we recommend to future workers on the model. The ability to calculate 
matrix elements of Sj between non-YHW states will be a key requirement. 



4 The particular application of our selection rule to matrix elements of the form (T| Sf |GS) where 
|GS) is the groundstate in a magnetic field, a — + and |T) is a YHWS, has now been confirmed 
mathematically in [0] using the theory of Jack-polynomials. Additional evidence for the Sf case can 
be found in those references as well. 
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3.6 Appendix. 
3.6.1 Appendix A. 

We show that localized spinon wavefunctions with N sp spinons all pointing up are 
necessarily linear combinations of YHWS with exactly N sp spinons. This follows 
easily from the fact that both S + and J± annihilate these states. The derivations 
is very similar to that of Appendix A of chapter |l[ If we write \l/(ni, . . . ,Um) — 
ip(z ni , . . . , z nM ) where z n = exp(^p) and 

M N sp 

ip(wi, . . .,w M ) = - WjfJlwi J[ Y[(wi - z a .), (3.23) 

i<j i i=l j=l 

then 

N N 
(S + ^)(n 1 ,...,n M ^ 1 ) = 53*(ni,...,7lAf_i,j) = V^m, ■ ■ ■ , Znu-u e ~^) 

3=1 i=i 
= <i[>(z ni ,...,z nM ,0) = 0, (3.24) 

(where the J2j was recognized as the zero- mode of a Fourier expansion). And with 

Wjk = —2——- ; 

3 K Zj-z k 

M N 

(J 1 + ^)(n i , . . . ,n M -i) w nj! i^(ni,. . . ,n M -i,i), (3.25) 

3=1 i+rij 

where we used the fact that 2Sf = 1 — 2 J2j=i Si, nj and that \I/ is a symmetric function 
which vanishes when two of its arguments coincide. Employing eq. ( 1.48 ), we have 



Ji^I = 0. since if> has a double zero when two of its arguments coincide and it 
vanishes at z = 0. 

3.6.2 Appendix B. 

In this Appendix we want to prove eq. (|3.11|) . Let us first introduce the following 
identity which holds for any set of distinct complex numbers {uJi}: 

M-l m m M 

n & - *) = e n ^hw n^* - <**)■ ( 3 - 26 ) 

i=l fc=i z(^fc) L ° k i=i 
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The RHS is just the Lagrange interpolation formula applied to the function in z on 
the LHS! 

Say we want to write Sf^io^x as a linear combination of localized spinon wave- 
functions with two more spinons than ^{ ai }, and all spinons pointing up. We fix one 
of the additional two spinons at i, the site on which acts, i.e. : 

(S+^jajXni, . . . ,n M -i) = $(ni, . . . ,n M -i,i\aci, ...,a Nsp ) 

= a p^( n i, ■ ■ ■ ,n M -i\®i, • • • ,a Nsp ,i,p). 

(3.27) 

Using eq. ( |3.8| ) we can divide out common factors of [z nk — z ni ) etc. , and we are left 
with: 

ZiY[(z nj - Zi) = a P Y[(z n] - z p ). (3.28) 

The result follows when we apply the identity (|3.26|) to this equation with z = Zi and 
a P — z i IIa: wnere ^ ne z k are randomly chosen distinct sites which don't coincide 
with the localized spinons. 

3.6.3 Appendix C. 

We show that Sf |GS) can only have or 2 more spinons than |GS) where |GS) is a 
YHWS groundstate in a nonzero magnetic field with 

(ni, . . . ,n M \ GS) = ^ (^i, • • • ,%) 

= n(^-^) 2 Il4" M+1 . (3.29) 

i<j i 

The proof hinges on the fact that in either of those cases (AN sp = or AN sp = +2) 
acting on 5"? |GS) twice with J± will annihilate that state. Potential AiV sp = 4, . . . 
contributions should survive as they are at least 2 levels from the top of their Yangian 
multiplet. Now 
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(2S?*o) 



M 

i=i 



(3.30) 



and 



2J+S?* ) (%,..., n M -i) 



M-l iV 

- E E «w x 

M-l 

-2 ^ u? nfc ,i*o(ni, • • • , n M -i, i 
fc=i 



M-l 



\E r o( n i.---,«M-i, 

i) - 2 E ^,n 3 ^o(niv..,n M -l,0 



(3.31) 



The first two terms vanish when we apply the convolution theorem with Wki, as in 
Appendix |3.6.1| ; only the last term survives. Notice that the {ni}; i = 1 . . . , M — 1 
cannot be equal to % anymore. This state is trivially annihilated by S + as it vanishes 
at z — 0. Furthermore: 



s M-2 V M-2 

Sf^oJ (ni,...,n M - 2 ) = EE w n P ,q E w n k ,i^o(ni, ...,n M -2, 

p=l q^=i,n p k=l 
M-2 V 

+ E E w n p , g ^g,i^o(^i, ■ • ■ ,n M -2,q,i)- (3.32) 
p=i g^«,n P 

With the help of the convolution theorem we can set the first term to zero (we can 
stick in the extra term with q — i, which is missing, at no cost since vanishes 
when q — i). For the second term we use the identity that lies at the heart of the 



integrability of the Haldane-Shastry type models: WijWjk + WjkWki + Wkitu 



-1. 



JZ) SfV) (m,...,n M -2) = 

M-2 N 

~EE (1 +Wi : n p Wn p ,q + W qi iW i:np )ty(n 1 ,...,n M ^ 2 ,q,i). 
p=l q^=i,n p 



(3.33) 



The first term is zero since \I/ is a polynomial that vanishes at the origin, whereas 
terms two and three can also be put to zero with the help of the convolution theorem. 
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Chapter 4 

New Types of Off-Diagonal Long 
Range Order in Quantum 
Spin-Chains 



4.1 Introduction 

In this final chapter we attempt to link properties of the ISE-model to some other 
one and two dimensional correlated systems. We want to show that there exists 
a new form of Off Diagonal Long Range Order (ODLRO) in ID spin models; this 
ODLRO, inspired by electron pair hopping in coupled Hubbard chains or planes, 
was suggested by Anderson in 1991 [0]. The corresponding order parameter also 
has an interesting analogy with the Read-Girvin-MacDonald order parameter in the 
Fractional Quantum Hall Effect (FQHE) [0, 0]. To be specific, consider the following 
linear operator A(j), acting between the 2 7V_2 -dimensional (N — 2)-site spin Hilbert 
space, and the 2^ dimensional iV-site Hilbert space. Acting at site j, it consists 
of slipping in 2 sites, with two spins in a singlet, between sites j and j + 1, in the 
groundstate of a particular model. This operator acquires a non-zero expectation 
value in the groundstate in the following sense [0]: 



(0,JV-2|At(j)A(0)|0,JV-2> 



J>, N - 2\&{j)\v, N)(u, N\A(Q)\Q, N - 2) 



= |<0, JVjA(0)|0, JV- 2> +... 
-> const foriV»j»l. 



(4.1) 
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|0, N—2) is the groundstate on N—2 sites, and {\v, N)} is a complete set of eigenstates 
of the N site model, where v — labels the groundstate. In the case of translationally 



invariant models, the first term in the expansion ( |4.1| ) guarantees the ODLRO for A(j) 
provided it exceeds ~; in that case the eigenstates have a definite momentum and 
(0,N -2\&(j)\v,N) = e^ p o~ 2 -P^(0,N - 2\& (0)\v, N) . The sum (0) is bounded 
from below by 2 (0, JV - 2|A(0)|0, N) 2 - 1. With the notation: 
a u = e i( -Po~ 2 -p^)j. |(z/, iV|A(0)|0, N — 2)| 2 , we find this, since \T v a v \ > |a |-£„>o Kl 
= 2 1 ao | — 1 if | ao I > \ an d \ a v\ = 1- The last equation is a consequence of the fact 
that A(0)|0, N-2) has norm 1. 

We tested for this form of ODLRO for three different ID models: the short ranged 
NNE and XY models, and the long ranged ISE-model. In all three cases we find 
the ODLRO for A, and (even more remarkably) the constant in eq. ( f4.1|) is almost 
insensitive to the specific model that we test for: ~ (0.81) 2 ~ 0.65. A(j) is mentioned 
implicitly in [0], but the context is that of a ID Hubbard model. We will turn our 



attention to that in section 4.3 



One might argue that in one dimension only quasi-long range order with algebraic 
decay of the order parameter correlation function should occur — true LRO being ruled 
out by the Mermin- Wagner theorem. The way out is that true LRO is allowed if the 
order parameter either commutes with the Hamiltonian, or if it breaks a discrete 
symmetry. The latter is the case here, we will find that the singlet insertion breaks 
the Z 2 symmetry of shifting by one lattice site. 

In the next section we will discuss numerical results and analytical calculations, 
indicating the presence of the new ODLRO. In section [4.2.2| , we will elucidate the 
analogy between the singlet insertion operator A and the Read order parameter in the 
v = \ FQHE. Finally, section |4.jj| will illustrate why the LRO for a singlet insertion 
may have important repercussions for superconducting correlations. 
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4.2 Evidence of Singlet Insertion ODLRO. 

4.2.1 Spin- and Spinon-Singlet Insertion in the ISE- model 
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Figure 4.1: Overlap between the (N — 2)-site groundstate with a singlet appended, 
with the iV-site groundstate, as a function of N, for three spin models. 

In fig. ( |4.1|) we have plotted numerical results for this singlet insertion 
(0, N\A(j)\0, N — 2) , for the groundstate of the three models mentioned earlier, 
with Hamiltonian: 



h = £ Ai - j) {\{stsj + s.sf) + a(s*s*)} . 

&3 



(4.2) 



For the XY and NNE models J(n) = 5 n> i + 5 n -i and A = 0, A = 1 respectively. For 
the ISE-model Jin) = , N . , nv ^ , A = 1. 

(^sHlT )) 2 ' 
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A^(j)A(0) can be rewritten as 

At( J )A(0)=P 1 5 2 (P 12 ..,) 2 (4.3) 
where Pf 2 is the projector for the spins on sites 1,2 onto a singlet: 4 — Sj • S 2 , and Pyi -j 



is the operator that cyclically permutes spins on sites 1 • • - j. With the help of eq. ( |4.3j ) 

2 



we have verified that 



<0,iV|A(0)|0,iV-2> 



is indeed the dominant contribution to 



(0|At(j)A(0)|0) for iV -> oo. 

These overlaps are shown as functions of A" in fig. [0] together with fits to the 
results of the form 0.817 + 0.778A^ 2 for the ISE model, 0.820 + .740A^ 2 for the NNE 
Heisenberg model, and 0.808 + 0.8 19iV -2 for the XY model. The results strongly 
suggest that all overlaps remain finite in the limit where the system size goes to 
infinity. Note that since the phases of the iV and N — 2 site wavefunctions may be 
chosen independently the phase of the overlaps is meaningless and further, since the 
groundstate momenta of the iV and N — 2 site groundstate wavefunctions differ by 
tt, the overlap is multiplied by minus one if the location of the singlet pair is shifted 
by one site. We have chosen the overlap real and positive for convenience. As a 
consequence of this alternation and the finiteness of the overlap, the singlet insertion- 
singlet deletion correlation function has fixed magnitude and a sign determined by 
the number of sites between the insertion and the deletion modulo 2 in the limit of 
infinite separation of the insertion and deletion; the models therefore have a hidden 
ODLRO which breaks a Z 2 symmetry. As fig. fl4.1| ) shows, the overlaps for the A^-site 
groundstate with the N — 2 site groundstate with a nearest neighbor singlet insertion 
in the ISE and Heisenberg models are not only finite in the N —>■ oo limit, but also 
surprisingly close to each other, despite the fact that the range of the interaction in 
both models is quite different. 

We now present an analytical calculation in the ISE model that gives some under- 
standing of origin of the finite overlap. The unnormalized groundstate wavefunction 
on A" sites of the ISE model in a basis labeled by the positions of the down-spins 
{ni, . . . , n M }, M = f is given by eq. (|Hj ): 



N N 
2 2 



(n l5 . . . ,n f ) = Ilsin 2 (^^tt) (4.4) 



i=l i<j 
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In a basis of local spins {|oi • • -a M )}, <Ji = ±| this reduces to [0]: 



M 



^>i--^m} = ii< 

i<j 



Zi — Zj) " i,a ^e 2 



(4.5) 



It: in 
N } 

-TTll'l 

> JV-2 



„ =1 , while for the 

N-2 
1 



For the N site ISE groundstate M = N and = Cn = {e 

iV - 2 site groundstate M = N - 2 and {z^} = CAr_2 = {e^}%=i- To compute 
the ISE overlap we add crjv-i, cat, sitting in a singlet, to the iV — 2 site groundstate. 
This overlap is not calculable since the Zj's from both sets are not commensurate 
with each other. However, if we slightly deform the set Cat_ 2 to be {e 2 ^}^ 2 and 

(2) 

leave (JN-i-> a N m a singlet then we obtain a new state W)y , that can be recognized 
as a localized 2-spinon state [0, 0], eq. (|3.8j ), as in fig. fl4.2|) . This is an uncontrolled 




Figure 4.2: Graphical illustration of the deformation of the singlet insertion to a 
spinon singlet insertion. The gray area is the groundstate wavefunction, or an adia- 
batic deformation thereof. 



approximation whose justification lies in the agreement with our numerical results. 
This state is an admixture of eigenstates that contain only or 2 "real" spinons. Here 
the localized spinons sit at sites N — 1 and iV in a singlet. In general we could have 
put them at sites a, p by deforming Gjv-2 m to {e n \% = 1, . . . , jy }/{e n , e n }. 
In the basis of states with M overturned spins with respect to the ferromagnetic 



state labeled by their positions along the chain: m, . . . , tim, anc l are given 



f (2) 
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by 1 : 



(n 1( 



nc-r n^m 2 



2 1 



JJ sin (^tt) sin l 



(4.6) 



i=l 



i=l 



The overlap between and \I/ jy* for arbitrary separation a — (3 between the spinons 
is calculated in Appendix A. We find that the normalized overlap is given by: 

2 sin (vr^) 



(4.7) 



AT sin(^) ^ Si(7r(a -/?))' 

where Si(x) is the Sine integral function. Thus for a nearest neighbor the overlap 
0.82917, which is within 1.5% of the Heisenberg and ISE groundstate 



is 2 



singlet insertion overlaps 



4.2.2 Relation to the FQHE. 

We noticed that the action of the generalized singlet insertion operation at sites a, /3, 
on the ISE-model groundstate was essentially the following: it creates two "holes" 
at those sites — where the down spins can't reside — by way of multiplying the wave- 
function by a factor Ili^n; — z a)( z m — Zp). Then it inserts a down spinon at site 
a, which, in the down spin - particle analogy is just a (hard-core) boson 2 . This all 
looks very similar to what we know about the bosonic v = | Fractional Quantum Hall 
Effect. In this theory, particles are constrained to a 2D plane in a very high magnetic 
field, such that the cyclotron frequency dominates over the interactions between the 
particles. In that case the Hilbert space is effectively just the set of lowest Landau 
level wavef unctions, of (kinetic) energy ^htu c = y— , which are degenerate 3 . When we 
turn on interactions, some dramatic things happen: when the filling fraction v = y 
is close to a rational number of the from — , the system develops a gap, and becomes 
an incompressible fluid with approximate groundstate wavefunction [0]: 



2 



...,z M ) = Y[( Zl - Zj ) m e-^, (4i 

i<j 



1 The conversion from {e^} to basis is based on the identity Oj=i 0- ~ T~) = cons t when 

= e^ri and z lies on the lattice z N = 1. 
Actually, we have the down spinon at a half the time, and (3 the other half, but just fixing it at 



a leaves the matrix element the same, up to a factor see eq. ( 4.22 ) 
3 Degeneracy — g — ■§- = -^x 2 ; A is the area of the sample 
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where I 2 = and Zi = Xi + iyi is the complex co-ordinate of particle i. Obviously, 
for m even the particles must be bosons, and fermions otherwise. Wavefunction (|4.8j ) 
is valid for particles confined to a circular droplet. 

Elementary excitations of this wavefunction are so-called quasiholes and quasipar- 
ticles. We will concentrate on the quasiholes. To create one at co-ordinate z, we have 
to introduce an extra zero in the groundstate wavefunction at that point: 

N 

(zi, ...,z N \z) = H(z - Zl ) tf (*i, • • • , z N ). (4.9) 

i=l 

This is mathematically identical to the up-spinon insertion in the ISE spin wavefunc- 
tion. A quasihole is a vortex in the quantum hall fluid, that allows one extra flux 
quantum to pierce the fluid. 

The entry into a Quantum Hall phase — with gap — can be signalled by the non- 
zero expectation value of an order parameter. The Quantum Hall order parameter 
is not of simple the Bose-condensation type ^(r), or Cooper pair type ipki^-ki- A 
first version was suggested by Girvin and MacDonald [0]. Soon afterwards, Read 
introduced his version [0]: 

6 m (z) = exp (m J d 2 z' \n(z - z')p(z')) ^(z), (4.10) 

whose effect is to introduce a neutral complex of m quasiholes, and a particle, at 

co-ordinate z. This operator has true long range order in a v = — Quantum Hall 

state. Thus, in the limit where the number of sites in the ISE-model becomes large, 

the spinon-insertion expectation value is essentially the Read order-parameter in the 

bosonic m = 2, v — | FQHE, acting at the edge of the circular droplet. 

The analogy can be extended further. When the spinon sites are dragged apart, 

such that 1 <C \a — /3| -C N, we are essentially measuring the spinon propagator. 

Recall that the localized two spinon wavefunction is a linear combination of eigen- 

states with two real spinons or none (the groundstate). When we overlap with the 

latter, we measure the amplitude for one of the spinons to propagate to the other and 

annihilate (provided that they are in a singlet configuration, since S 2 is conserved). 

The decay of this propagator (0, N\(Q, a); N — 2) (where the ket designates the two 

_i 

spinon state) can be derived from eq. (|4.7| ). It is proportional to a 2; l<a<JV. 
We can do the same thing for the singlet insertion, by not inserting the spins in the 
singlet next to each other, but a distance a apart. This correlator is plotted in fig. 
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( fO) , and shows the same decay oc a 2 . We have evidence for the same exponent in 
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Figure 4.3: The overlaps between the N site groundstate and the N — 2 site ground- 
state with two spins in a singlet inserted at sites 1 and y. XY model results are from 
exact numerical results, whereas the ISE-model data for iV > 18 is from Monte Carlo 
simulations. 



the NNE model 4 . 



4 Actually, fig. (4.3) shows the decay as a function for the system-size N, not a. The reason is 
that for periodic boundary conditions, if the matrix element decays algebraically, we expect it to be 
of the form 



(0,N\(a,0);N-2) cx 



1 



^sin| 



7TCX ^ 

N , 



2 A 



(4.11) 



That is, we assume that the ISE-modcl is conformally invariant at T — and the correlators of 
the primary operators in this theory, including the singlet insertion, decay algebraically on the 
infinite (1+1) dimensional space time: (<p(z)(j)(0))) ~ M~ 25 , and z = x + it. With the mapping 
w = u + iv = y~ ln(z) the spatial direction v becomes compact, and at equal times (u — 0) we obtain 
the desired result [0]. We will see later that A(r) is indeed a primary operator in a c = 1 CFT. 
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For the v = | FQHE, the corresponding quasihole propagator decays exponentially 
in the bulk of the material (because of the gap). Low energy excitations can propagate 
along the edge of the fluid, which is contained by a confining potential. Such modes are 
called edge states. Quasiholes at the edge have a propagator that decays algebraically, 
and with Kac-Moody algebra techniques, Wen [0] showed that its exponent is | as 
well. As a remarkable aside: the exclusion statistics of the spinons [0], which is | 
(bosons: 0, fermions 7r), is the same as the exchange statistics of the quasiholes in 
the v = \ effect [0]. 

4.2.3 ODLRO in the XY-model. 

In an effort to better understand our results for the Heisenberg and ISE models, we 
examined the overlap for the same insertion operator acting on the N — 2 site XY 
model groundstate with the N site XY model groundstate. The strikingly similar 
result, here extended to a larger system size than was possible for the Heisenberg and 
ISE models, is shown in fig. |4.1| . Since the XY model can be mapped onto spinless 
fermions, the "order parameter" can be recast in that language and a more detailed 
study made. In that language the non-zero overlap is between the groundstate for 
N/2 spinless fermions, with either periodic of antiperiodic boundary conditions, on 
N sites (with positive hopping integral so that their momentum in the groundstate 
are centered around tt) and the state obtained by adding two sites and one fermion 
(in a superposition of being on the two added sites with a relative minus sign for 
the two different sites) to the groundstate for the y — 1 spinless fermions on iV — 2 
sites with the opposite boundary conditions from the N site case. The change of 
boundary conditions is a non-local operation in the fermion language, however in the 
spin language it arises from a local operation, i.e. the insertion of two additional spins 
in a singlet configuration. 

Since the wavefunctions for the two fermion states are just Slater determinants, we 
can bound the singlet insertion matrix element rigorously from below. The technical 
details can be found in Appendix B. It is found there that \(0, N\A(j)\0, N — 2)\ 
> e~2 for N — > oo. 

It is somewhat remarkable that the ODLRO for the XY-model is the result of 
three conspiring effects, every one of which by itself would lead to a vanishing matrix 
element as iV —>■ oo. 
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• |0, N) had y electrons, all in extended states and A(J)\0,N — 2) has y — 1 
electrons in extended states too, but the last is in a localized two site state, 
resulting in a decay oc 

• The extended states in the iV-particle groundstate have periodicity N, whereas 
the (N — 2)-particle groundstate has periodicity N — 2, with a resulting orthog- 
onality type catastrophe in the overlap. 

• The particles in the iV-site groundstate have periodic/antiperiodic boundary 
conditions, whereas the (iV — 2)-site version has antiperiodic/periodic ones (de- 
pending on the parity of y). 

For the XY-model too, we have done numerical calculations, studying the decay of 
the singlet insertion as a function of the separation of the singlet spins, see fig. ( |4.3|) . 
In this case we find an exponent very close to |. In the next section we will give 
arguments for this particular decay exponent from the point of view of the Luttinger 
liquid picture, which describes the low energy excitations of these spin models [0] . 

Finally, the XY model and the ISE model both have wavefunctions of the Jastrow 
product form of the Calogero-Sutherland model [0, 0]: 



N JV 

(m, . . . , n f ) = f[(-r II sin- f^%) (4.12) 

i=l i<j v / 

For the XY model p — 1 whereas for the ISE model p = 2. We have studied nu- 
merically the overlap for a number of values of p and find that it remains finite in 
the limit of large system size and is essentially constant for 1 < p < 2. We have 
also examined < p < 1 and find the overlap to be finite in the large system size 
limit and monotonically decreasing with p down to a value of 1/2 at the point p = 
where the wavefunction is that of Gutzwiller projected free bosons and the can be 



computed exactly as 77^ y j^zf- Note that this implies 0(-^)) corrections in contrast 
to the O(-jkn) corrections we find for the XY, ISE and Heisenberg models. It is not 
clear when corrections of size O(jt) first appear. 

It may be tempting, based on the similarity of the wavefunctions Q4.12 ) to relate 



the XY model ODLRO to a v — 1 Read order parameter, acting on the edge. However, 
the singlet decay exponent | is not consistent with the electron propagator exponent 
in the Quantum Hall situation. 
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4.3 Singlet Insertion in the ID Hubbard Model. 

In the context of spin models, the singlet insertion operation may seem somewhat 
esoteric, requiring the Hilbert space to change size. Now if we imagine the sites to be 
physical electrons, then the singlet insertion procedure corresponds to hopping two 
electrons in a singlet configuration onto a chain. 

We follow the tradition of studying spin models in the guise of the ID Hubbard 
model in the limit of one particle/site, and U — ► oo [0], and consider the action of our 
order parameter in the context of this simplest of interacting ID itinerant electron 
models. Its Hamiltonian is given by: 

# = t5>Lc i+lff + h.c.) + C/£> iT n a . (4.13) 

is i 

The i's are site-labels. The U > term discourages electrons to occupy the same site. 
There is an additional non-academic reason to study the effects of singlet insertion 
in a Hubbard-type model. One of the mechanisms to explain the high transition 
temperatures in the copper-oxide superconductors is built around the relevance of 
the hopping of electron singlet pairs between the 2D layers in these materials [0]. 
We will attempt to decide on the relevance of this operator in the one-dimensional 
Hubbard model for U — *■ oo, based on the ODLRO of our singlet insertion (see [0] for 
possible ways to extend this procedure to two dimensions). The inter- chain hopping 
operator acting at sites j, j + 1 is given by: 

= O l (j)Ot(j) +h.c, (4.14) 

where the electron creation operators for chain 1, and the gTs act on chain 

2. The relevance of this operator in the renormalization group sense for perturbing 
away form the uncoupled chain case, is determined by its scaling dimension d s in 
~ j 2ds for j -> oo. If d s > 2, singlet pair hopping is irrelevant. |0) 
is the tensor product groundstate of the two uncoupled Hubbard models. The Os-Os 
correlator then factorizes into an (O\(j)Oi(0)) and (.7)^2(0)) part. Ogata and 
Shiba have pointed out an additional simplification for large U [0]: the groundstate 
of ( |4.13| ), for M electrons, factorizes into the convolution of a spin and a charge part: 



^0(^1, • • • , n M \ai, a M ) = ^sf(«i, • • • , n M )x{ cr i, vm), 



(4.15) 
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where \1/sf is the determinantal groundstate wavefunction of M free spinless fermions, 
labeled by co-ordinates x(°i> ■ ■ ■ > °m) is the groundstate of the NNE Heisenberg 
model (at least for = 4m + 2, m integer) of M spins, residing on the spinless 
fermions. Then (0i|(9}(j)Ci(0)|0i) reduces further to: 

(*SF|^-+iV'iV'M(-) E<W)rH |*SF> • |(x|A(j) + A(0)|x>| (4.16) 

A is the singlet insertion order parameter from the spin models which alternates with 
j, i.e. keeps track of the number of fermions to the left of the insertion site, whence the 
(_)S,(<j) n! f ac tor in the spinless fermion matrix element. ipi,4>i are spinless fermion 
operators. The second factor in eq. (|4.16| ) is constant (~ 0.65 as we saw before). 
The matrix element decay eq. ( 4.16Q then follows from the fact that free fermion 



operators have dimension 1 and the alternating term has dimension ^. We want to 
derive this a bit more extensively in a bosonized language in the original Hubbard 
model. We will also include a spin anisotropy term to accommodate the XY model 
limit later on. In the continuum limit the Cj S — > ipR tS (ja)+ipL tS (ja), the field operators 
acting at the left and right Fermi points, s = ±^ is the spin label, and a the lattice 
spacing. In one dimension they can be rewritten in terms of boson fields [0, 0]: 



L "-0^ 

= **"-™«&\ (4 . 17) 

a is a cut-off scale. The fields (p^ and d x 6 P) p = p,cr are canonically conjugate bosonic 
fields: [(p^(x), d x 9 u (x')] = m5^ u 5{x — x'). One also writes ILj = 7rd x 9^. The fields with 
label a create/annihilate spin excitations (destroy up- and create down-spin electrons 
and vice- versa); the p- fields relate to charged excitations. A peculiarity of the ID 
Hubbard model is that these fields propagate independently: spin-charge separation, 
as is evident from the Hamiltonian in the bosonic language: 



2 

7TV Q K, 



h = i dx\^ui+^(d x( p P r 



+ Jdx {^f^K + ^(d x ^) 2 +$cos(>/8^)} + . . . (4.18) 

The ellipses refer to other irrelevant terms originating from the non-linear electron 
dispersion relation and the underlying lattice. v c and v s are the velocities of spin and 
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charge excitations respectively. K p and K a are functions of the Hubbard parameters. 
All correlation function exponents can be expressed in terms of these two. The spin 
cosine term is irrelevant if < U < oo and the anisotropy A satisfies — 1 < A < 
1. The non-linear term only acts to renormalize K p and K a : K p (U = 0) = 1, 
K P (U — > oo)= |, and K a (A, U — > oo) = 1 — - arccos(A) [0]. So at long wavelengths, 
we have two uncoupled free linearly dispersing bosons with different velocities. The 
electronic correlation functions can be deduced from the bosonic ones via ( |4.17| ). 
The latter are easy to compute, since the fixed point Hamiltonian is quadratic. For 
example [0]: 

\ a/3 



(JaORn {x,t) e -if3e Rll (0,0) \ 



Km 

4 



(4.19) 



Xx - v) 4 ^ 2 ( x + V) 4 M ; 

Notice that for interacting electrons (K p ^ 1) a right Fermi point electron operator 
has an anMiolomorphic part too! 

In the following we will neglect the phase-factors, since we are just interested 
in the scaling dimensions. Consider the case where the singlet electrons are in- 
serted, not as nearest neighbors, but at sites x, x', and extracted at sites y,y': 
(0\if}^(x)ip^(x')tp^(y)ijji(y') |0) + . . ., where 1 <C \x — x'\, \y — y'\ <C \x — y\ in units 
of the lattice spacing a. In this limit we can then apply the operator product ex- 
pansion of fields ij^(x) and ^(x') and separately ip(y) and ip(y'). In terms of the 
bosons: 



e ie m (xt) e ie Rl (x't') 



— 

e V2 



(Ax - v c At) Pp+ ^(Ax + v c At) Pp 



(xt) 



r(x't') 



| . e iV28 Rp (xt) . 

-:1 + ...:, (4.20) 



(Ax - v s At) Pa+ i (Ax + v s At) 
where p p = -^r- + and p a = + We used the generalization of the Glauber 

formula: e A e B =: e A + B e (°\ AB + A % B l°> for normal ordering with respect to the bosonic 
vacuum, and A, B commuting with their commutators [0]. There is a term with 
identical exponents for R — > L, and less relevant terms, containing mixtures of both. 
Notice that the equal-time singlet insertion decay exponent of fig. ( |4.3| ) can be read 
off from QO01) : 2p a = ^ + = \ for A = 1 (NNE-model limit) and § for A = 
(XY-model limit). 
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The part of the electron singlet operator that yields the spin-singlet correlator is 
e V2 D «- w . This operator is semionic, in that commuting two of these through each 
other results in a phase , for the isotropic Hubbard model. This is compatible with 
the earlier observation that the singlet insertion is essentially a two-spinon creation 
operator. [0] suggests a semionic operator in the context of a CFT realization of the 
low energy part of the ISE-model, bilinears of whose Fourier modes create 2 spinons 
as well. An identification of the boson-exponential with this CFT operator seems 
appropriate. 

From eq. ( |4.2(J| ) we learn that the electron singlet insertion in bosonized form is 
e iV26 Rp ^ w ith dimension | and correlator decay oc (x — y)~z, as we found before (see 
also [0, 0, 0]). The expression for O contains a charge-prefactor which grows with 
separation between the singlet spins, since it has a negative scaling dimension. This 
need not concern us, since this separation is bounded from above by the distance to 
the singlet extraction, i.e., the OPE is no longer valid when x' ~ y,y'- 

The decay exponent | is insufficient to make singlet pair hopping between chains 
relevant. From our spin model singlet insertion, we find the culprit to be the operator 
(— )^(<j) n \ caused by the alternation of the spin order parameter. It has been 
suggested [0] that this can be remedied if it could be canceled against a similar 
alternating factor in the second chain. 

4.4 Conclusion. 

We have shown that the operator that inserts a singlet into a ID spin chain acquires 
LRO in the ISE-, XY- and NNE-models. We provided numerical evidence and sug- 
gested that its LRO might be expected based on its interpretation as a two spinon 
insertion, as well as via an analogy to Read's order parameter in the v — | QHE 
(for NNE and ISE- models). The fact that the order parameter breaks a Zi symmetry 
causes the more physical electron singlet hopping operator in coupled ID Hubbard 
chains to become irrelevant. Bosonizing the anisotropic Hubbard chain allows us to 
calculate the decay exponent of the spin singlet insertion operator or, equivalently, 
the spinon propagator, as a function of the anisotropy. 

The fact that this ODLRO is present in three quite different antiferromagnetic 
isotropic spin-chains, raises the question whether there are models where it doesn't 
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occur. XY-, ISE- and NNE- models are all gapless and one might imagine that that is 
a key feature to this particular ODLRO, but for instance the Majumdar- Ghosh model 
with a gap above its dimerized groundstate possesses the singlet-insertion ODLRO 
by definition. This is an important issue since the usefulness of the order-parameter 
depends on its ability to discriminate between different models. 
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4.5 Appendix 

4.5.1 Appendix A: The Two-Spinon Matrix Element 

We compute the (normalized) overlap between the iV-site ISE groundstate \^/^) and 
the iV-site two-spinon wavefunction, localized at sites a and 0, \ty( 2 \a)). The two- 
spinon wavefunction is an adiabatic deformation of the (N — 2)-site groundstate plus 
a singlet. Because of translational invariance we can fix f3 to be 0. What we need to 
compute is: 



= A-B~\ 



(4.21) 



(fr(Q)|fr(2)( a )) _ (fr(°)|fr( 2 )(q)) / (y( 2 \a)\¥ 2 \a)) ^ 

^/<*(0)|^(0))^(2)( a )|^(2) (Q)> " (*(0)|*(0)) { (^(0)^(0)) J 

First we will determine (^ (0) |^ 2 (a)). Setting M = f - 1: 

(*(°)|*( 2 )(a)> = (0 £ 4* (2) K...,n M |0,a)x 

ni,...,n M 

{*^(m, . . . , n M , 0) - ^ (0) (n l5 . . . , n M , a)}* 
= )! a(a) (1 -(-)») . (4.22) 



Here 



e n 

ni,...,n M \i<j 



<T((l) = > | | | S11T f ~ ~ 7T 



M /vrnA il . /n ? -a 



x 

1=1 v / j =1 



n-'^j n-iv 1 ")' (4 - 23) 



In eq. (4.22) we shifted all the n; by a in the second term to bring both terms in the 
same form. 

Since is a singlet, we know that if the two spinons in were in the 5*2 = 
triplet state the overlap should be zero, i.e. (^r(°)|^( 2 ) (at), triplet) oc a(a) (1 + (— ) a ) = 
for all a. Therefore a(a) vanishes for all even a. At the same time we see from eq. 



( |4.22f ) that <j(a) is a polynomial in cos Ojf) • This is easily checked by expanding the 



sin { ni N a T^j and noting that a(a) is even in a. We conclude immediately [0]: 

JV 

a(a) oc n(cos(^)-cos(^)) 
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sin ( 


' not \ 
V 2 J 


sin | 


' not \ 


\ N ) 



oc y^-. (4.24) 

sin (t) 

which holds, even when a is not integer. This gives us the a dependence in A. To 
get the normalization correctly, let us write out A (a odd): 



.4 



1 (?) 


|l • 2 • V 

' • t—tn\ ,...,n 


M n 4<i sm 4 ( \ j tt) 


iniSin 3 (^7r) sin(^Tr) 


V2 


(?) 


l!-En 1 ,...,n M+i ni<jSm 4 ( \ 3 11 j 


l 



2 Em B]tf sin 4 figyr) n, sm 3 (g tt) sin to 



7T 



V2iVE ni> ..., nM Ui<j sin 4 (^tt) n, sin 4 ( 



N 



(4.25) 



In the last step we changed the summation variables rii — > +um+i] i — 1, . . . , M in 
the denominator and did the trivial um+i summation. Despite the fact that physically 
a is actually an integer, we can formally take the limit a — > in the last expression 
and obtain: A (a = 0) = Therefore: 



Now we turn our attention to the second piece. Only its numerator is a dependent: 

(^ (2) («)|^ (2) («)) = 

' N --l)\ £ h i2 \n 1 ,...,n M \0,a) n 



{ni,...,n M } 
v z ' m,...,n M i<j v iv 

n^(^)sin=(^. (4.27, 



But this sum can be recognized as the (S z (a)S z (0)) static correlation function in the 
ISE model [0], or — after an (exact) conversion of the sums to integrals (see chap- 
ter ^) — as the one-particle density matrix in the Calogero-Sutherland model at half 
filling [0, 0]. The result for large N is proportional to Sl j^ , where Si(x) is the 
sine-integral function Jq dx' sin ^f ' . To normalize B properly we write it out: 



B 



(?-l) 




„„ni<j sin 4 It sin 2 (>) 


i sin 2 (^vr) 


1 


(?) 


l! •En 1 ,...,n M+i n i < j Sin 4 ( \ >7t) 


1 
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2 n <<3 - sin 4 (ggj II sin 2 (g vr) sin 2 (sgj 

iV 2 E m ,...,n M n,<, sin 4 ft sin 4 

In the last step we shifted the summation variables, as we did previously for A, to 
simplify the denominator. We now recognize that for a = 0, B = 2/N 2 . As for A, 
the normalization can be obtained from the case a = and we find: 

B - wi£< (429 » 

l\ z ira 

The entire expression for the overlap, with normalization, then becomes: 

l^M^ (4 30) 

N sin(f ) VSi(Tra)' 1 } 

4.5.2 Appendix B: Lower Bound on the XY ODLRO 

We want to show that (0, iV|A(j)|0, N—2) is finite when iV — > oo. We begin, following 
Anderson's treatment of the orthogonality catastrophe [0], by writing the overlap as 
the yXy determinant of of the matrix formed from the inner products of the single 
particle wavefunctions comprising the two states: 

overlap = \M id \ = |(^ 1} |*J 2) )| (4.31) 

where, for y an odd integer and periodic boundary conditions for the N site case: 

iV+2(l-o) / x„ 

|*«> = £ , > ) (4.32) 

n=l y/N + 2(l-a) 

except that 



= -L(\N - 1) - \N)) (4.33) 



2 



and 



The z's are given by: 



*%\N-1) = 

2 

(*g f |iV) = (4.34) 



,2fi,, iV 1., 
Zfc,i = ex P(^v"( fc + "4 ~ 2^ 

. 27TZ iV „ ss / . 

^,2 = exp(-^— + (4.35) 
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The entries in the matrix are then given for j ^ y by: 

1 + z 2 

Mi j = z j>2 1,1 (4.36) 

y /N(N-2)(z i , 1 -z jfl ) 

and by: 

M t4 = ^ • (*,!)"". (4.37) 
The cofactors associated to the elements M- n are determinants of matrices of the 

l > 2 

form Xij = / l % , where fi and Aj, depend only on i, while gj and Bj depend only 
on j. Such determinants can be rewritten using: 

1*1 = mu ^^'f^;^ (4.38) 

Neglecting phases multiplying it, this leads to a rewriting of the whole determinant 
as: 

l M l = ITi(l + 4i) fafei ~ gy^j ~ z h2 )] 

\/2N[N(N -2)]— 11^(^,1-^,2) 

\ - ^m,l - 1 Ili(^m,l - Zi t2 ) , 
X _2 , 1 n _ r ^ ^-Oyj 

We are unable to compute this exactly, but we have shown by careful consideration 
of the infinite products that to leading order in jj the determinant is bounded below 

5 

by e~2 and must therefore remain finite in the limit of N — > 00. To illustrate, we 
examine a single term in the sum over m, T m , defined by: 

JV-2 

rp _ Z rn,l — 1 IIi= 2 l (Zm,l ~ Zj j2 ) (AA(\\ 

Z ^ + 1 H7 =Wm (^,i - *,i) 

using the fact that the z^\ are roots of unity and z ij2 of (-1) this can be rewritten as: 

rp _ Z rn,l ~ 1 T[j{z m ,l + Zj,l) (4 41) 

N Z m \ 

Removing the term in the numerator given by z m ^ — 1 and relabeling the remaining 
Zi t \ so that they are in one to one correspondence with the Zj t2 and j runs from 1 to 
^ we get 

T m = ^ " 1)2 n — i — (4-42) 
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Since both sets of z's come in complex conjugate pairs or are real, for every m there 
is an n such that T m = T*. Therefore, only the real part of T m matters and, since the 
product can be shown to be real, we may replace by 2. The prefactor of the 

product is then of uniform sign and the sum can be bounded by replacing the product 
with a lower bound. The smallest individual term in the product can be calculated 
and is given by 1 — -| + 0{-^) so that the product is bounded below by - as is T m . 
The prefactor of the sum over m is also real and bounded below by by similar 
arguments so that the whole determinant is real and bounded below by e~i. 
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Conclusion 



We have seen in chapters [l] and |2| that it is possible to show full integrability of the 
ISE-model, both in the sense of the existence of a transfer matrix that satisfies the 
Yang-Baxter fundamental commutation relations (just as for all other integrable spin 
models) and of the presence of a set of constants of the motion, constructed by taking 
the static limit of a dynamical deformation of the spin model. We found that to 
every eigenvalue of the ISE-Hamiltonian belongs a set of degenerate wavefunctions, 
that derive from acting on a polynomial-type Yangian Highest Weight State with the 
Yangian algebra. The Yangian transfer matrix representations told us unambiguously 
the spin content of this Yangian representation, that is, the state counting rules of 
the model's quasiparticles: the spinons. Those counting rules are non-trivial, leading 
to the fractional statistics nature of those particles, a phenomenon presaged by the 
fact that for two Yangian representations Ti and T 2 , Ti®T 2 9^2® Ti, in general 
(lack of cocomutativity). 

As far as the integrals of motion are concerned, we saw that such a set, that con- 
tains the Hamiltonian i^isE could indeed be constructed, but this required a somewhat 
unnatural deviation through the dynamical deformation of the static spin model. A 
more direct approach should be possible once we can solve the model involving the 
introduction of a second quantum parameter, deforming it in the direction of the 
hyperbolic models, i.e. introducing an exchange integral in i?isE which is an elliptic 
Weierstrass function P(x\uji, 0J2) [0, 0]. Since the NNE-model is obtained in the limit 
uji = iN] u 2 — > 0, and the ISE-model when ui = N; u>2 — > 00, we could interpolate 
between the traditional Nearest Neighbor constants of motion to the more esoteric 
ISE ones. 

In chapter || we stated the remarkably simple, empirically found, selection rules 



106 



Conclusion 



107 



for excitations that can be generated by action of a local spin operator. The rule 
essentially states that the only Yangian multiplets that are connected via Sf, a = 
1, 2, 3 are essentially those whose motifs differ little: the number of their Drinfel'd 
zeroes in any range [a, b], 1 < a < b < N differs by less than two. This rule has 
only been proved in the special case of YHWS in both ket and bra of the Sj matrix 
element, via Jack polynomial Clebsch-Gordan coefficients [0]. The question of the 
extension to the case of non-YHWS matrix elements is still open. 

As a corollary to the Yangian selection rule, Sf acting on the groundstate in a 
given magnetic field, changes the number of spinons by no more than two, and can 
add up to two magnons (which move in opposite directions). In low fields, the dom- 
inant contribution to the spectral function comes from the addition of two spinons, 
in high fields, from the addition of one magnon. These selection rules predict the re- 
gions in the energy-momentum diagram that can be accessed in a neutron scattering 
experiment. The actual intensity of the inelastic scattering cannot be determined for 
all possible excitations. Most notably, (less polarized) excitations further down in a 
Yangian multiplet are not available. One would have to know what the permutation 
algebra of Yangian lowering operators with local spin operators would be, to be able 
to compute a typical matrix element like (r|T + (Ao)S|0). If one believes that contri- 
butions from states further down in a Yangian multiplet are related to those involving 
their parent YHWS |T) — as is the case for excitations of types (i), (iv) and (v) b in 
table |3J] — it seems reasonable to insert diagonal elements of the transfer matrix be- 
tween Si and T+ in (r|T + (Ao) T °A (u) Si|0). Since |0) is known analytically, one could 
choose T aa judiciously, such that its action on Sj|0) is manageable. Commuting it to 
the left through T + (Ao) will fix the matrix element. 

Chapter ^ presented evidence for the surprising existence of a new from of ODLRO 
(off diagonal long range order) in the groundstate of one dimensional antiferromag- 
netic spin chains, that consists of inserting a singlet into the groundstate, and ex- 
tracting it somewhere else. The ODLRO is present in short ranged XY and NNE and 
long ranged ISE spin models, with virtually identical values for the order parameter, 
which alternates in sign depending on where in the chain it is applied. This raises 
the question which models would not have this ODLRO. The ODLRO could be seen 
to exist, at least in the ISE case, by the fact that it is related to the spinon propa- 
gator, which can be shown analytically to have a finite piece in the thermodynamic 



Conclusion 



108 



limit. The spinon propagator bears resemblance to the Read order parameter in the 
v — \ Hall effect, which is known to have true long range order. The singlet insertion 
operation was applied to coupled one dimensional Hubbard chains, with large on-site 
repulsion. The corresponding operator makes a singlet pair of electrons hop from one 
chain to the other. That operator does not have LRO, in fact it is not even relevant 
(its correlator decays too fast), because the spin model order parameter breaks the Z 2 
symmetry of translating by a lattice site. Relevancy could be restored if the chains 
are properly correlated, such that the symmetry breaking alternation term in both 
chains cancels. 

The big question here is if, and how this can be extended to two dimensions. 
Numerically, things get difficult very quickly as a function of the system size, and 
there is no two dimensional ISE model that we may borrow intuition from. The issue 
of how to insert the singlet in a natural way into a 2D lattice is unclear as well. One 
could also consider generalizations to models with higher spin and SU (n) spins. 
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